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Abstract: A new model (CUJET3.0) of jet quenching in nuclear collisions coupled to bulk 
data constrained (VISH2+1D) viscous hydrodynamic backgrounds is constructed by gen¬ 
eralizing the perturbative QCD based (CUJET2.0) model to include two complementary 
non-perturbative chromodynamical features of the QCD confinement cross-over phase tran¬ 
sition near T c ~ 160 MeV: (1) the suppression of quark and gluon chromo-electric-charged 
(cec) degrees of freedom and (2) the emergence of chromo-magnetic-monopole (cmm) de¬ 
grees of freedom. Such a semi Quark Gluon Monopole Plasma (sQGMP) microscopic sce¬ 
nario is tested by comparing predictions of the leading hadron nuclear modification factors, 
R\ A (PT > lOGeV/c, y/s), and their azimuthal elliptic asymmetry vlf {pr > lOGeV/c, yfs) 
with available data on h = tt, D, B jet fragments from nuclear collisions at RHIC(-/s = 0.2 
ATeV) and LHC(y / s=2.76 ATeV). The cmm degrees of freedom in the sQGMP model 
near T c are shown to solve robustly the long standing Raa vs V 2 puzzle by predicting a 
maximum of the jet quenching parameter field q(E,T)/T 3 near T c . The robustness of 
CUJET3.0 model to a number of theoretical uncertainties is critically tested. Moreover 
the consistency of jet quenching with observed bulk perfect fluidity is demonstrated by 
extrapolating the sQGMP q down to thermal energy E ~ 3 T scales and showing that 
the sQGMP shear viscosity to entropy density ratio 77 /s ~ T 5 /q falls close to the unitarity 
bound, 1/47T, in the range (1 —2)T C . Detailed comparisons of the CUJET2.0 and CUJET3.0 
models reveal the fact that remarkably different q(T) dependence could be consistent with 
the same Raa data and could only be distinguished by anisotropy observables. These 
findings demonstrate clearly the inadequacy of focusing on the jet path averaged quantity 
(q) as the only relevant medium property to characterize jet quenching, and point to the 
crucial roles of other essential factors beyond just the (q), such as the chrorno electric and 
magnetic composition of the plasma, the screening masses and the running couplings at 
multiple scales which all strongly influence jet energy loss. 
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1 Introduction 

In ultrarelativistic heavy-ion collisions at the BNL Relativistic Heavy Ion Collider (RHIC) 
[1-4] and the CERN Large Hadron Collider (LHC) [5-12], strongly-coupled quark-gluon 
plasmas (sQGP) are created with deconfined color degrees of freedom under extremely hot 
conditions at sufficiently high beam energies [13-15]. This new phase of Quantum Chro- 
modynamics (QCD) matter provides a unique environment that contains an abundance 
of information about both the perturbative and nonperturbative aspects of QCD. In such 
collisions, large transverse momentum (pr) partons are produced in the pre-thermal stage 
and subsequently traverse the entire medium. They scatter strongly with the dynamical 
degrees of freedom in the thermal medium, and undergo both radiative and collisional en¬ 
ergy loss [16-22]. This “jet quenching” effect leads to suppressed yields of high px hadrons 
in nuclei-nuclei (A+A) collisions compared with the yields in scaled proton-proton (p+p) 
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collisions at the same center of mass energy. It has been discovered experimentally at 
both RHIC and LHC and serves as one of the key evidences for the formation of sQGP 
in heavy-ion collisions [13-15]. Jet quenching observables probe the color composition and 
chrorno field fluctuations in high density QCD matter that can provide insight into the 
novel dynamical mechanisms responsible for the surprising “perfect fluidity” property dis¬ 
covered through the systematics of bulk collective low px azimuthal harmonics at RHIC 
and LHC. 

Among a variety of jet quenching observables, two primary informative ones (at single 
hadron level) 1 are the nuclear modification factor Raa which is defined as the ratio of the 
A+A spectrum to the p+p spectrum scaled by the number of binary collisions, 


Raa(pt,v;Vs, b ) 


dN aa/ dyp T dp T 
Nbin dN p p/dyp T dp T ’ 


( 1 . 1 ) 


and the azimuthal elliptical harmonics V 2 which is the second Fourier coefficient in 


dN 

dypxdpT dcj) 


( PT,(/>,y;Vs , b ) 


1 dN 

2ir dypxdpT 


OO 

1 + 2 v n cos (n(0 - VBQ)) 

n =1 


( 1 . 2 ) 


They characterize the overall magnitude and the angular dependence of jet quenching in 
heavy-ion collisions, imposing stringent constraints on the mechanism of parton-medium 
interactions in jet energy loss models. However, most perturbative QCD (pQCD) based 
jet quenching frameworks have been found to fail to describe simultaneously 2 the high px 
light hadrons’ and open heavy flavors’ Raa and V2 data at RHIC and LHC [24-28]. 

One of the main problems with the conventional perturbative QCD picture of the 
(HTL) quasiparticle degrees of freedom in Quark Gluon Plasmas is that leading order 
(LO) pQCD estimates of the QGP’s shear viscosity to entropy density ratio [29-32], 

1 „ °- 071 . (1.3) 

s a 2 s \og{l/a s ) ’ 

predict that this ratio should be of the order unity, which is an order of magnitude larger 
than the unitarity p/s = 1/4 tt lower bound [30, 33] that was found to be required to ex¬ 
plain the “perfect fluidity” property of QGP produced in high energy nuclear collisions [13, 
14, 25, 34-36]. While the factor of ~ 5 quenching of hard leading hadrons observed in 
central collisions with Raa ~ 0.2, was well predicted [17] even with perturbative QCD jet 
medium coupling, the collective bulk azimuthal flow moments observed at RHIC and LHC 
appear to require much stronger interactions such as those assumed, e.g., in AdS/CFT 

1 In the present study we focus on single hadron observables, and the current CUJET implementation 
considers energy loss of single partons that are subsequently mapped to hadrons. This should be dis¬ 
tinguished from studies that focus on full jet evolution for describing reconstructed jet observables. The 
extension of CUJET framework toward full jet study would be a future project. 

2 It should be pointed out that recent studies on full jet observables (see e.g. [23]) provide quite reasonable 
descriptions of reconstructed jet Raa and V2 measurements at the LHC. The different situation for single 
hadron versus full jet results may be due to different sensitivities of different types of observables, or may 
also plausibly hint at systematic uncertainty in various studies of single hadron observables due to limited 
understanding of hadronization processes. 
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black hole modeling of the sQGP to account for perfect fluidity [25, 33] The long-standing 
“jet” Raa{pt > 5 GeV) vs “bulk” V2(pt < 2 GeV) [24, 25] as well as the “jet” azimuthal 
V2(pt > 10 GeV) puzzles [27, 28, 37, 38] continue to raise critical questions about our un¬ 
derstanding of the dynamics and composition of the QGP medium produced in high energy 
nuclear collisions and especially the consistency of information derived from high px > 10 
GeV jet quenching observables and low px < 2 GeV bulk collective flow observables at 
RHIC and LHC. While the soft hadrons’ V 2 (originating from collective flow of the nearly 
perfect bulk fluid) and the hard hadrons’ V 2 (due to geometric anisotropy of in-medium 
path lengths) are phenomenolgocially via different mechanism, they are both generated 
by the same underlying QGP medium for which a consistent microscopic model should 
account for both the soft and the hard transport properties. Many critical questions need 
to be addressed here. How do the effective degrees of freedom of nonperturbative QCD 
origin blend interpolate between the confined Hadron Resonance Gas (HRG) world at low 
energy energy density into an asymptotically free quark gluon plasma at extreme densities? 
Are there effective quasi-particles in the highly non-perturbative non-conformal tempera¬ 
ture range near the deconfinement transition temperature T c ~ 160 MeV? How do those 
“emergent” degrees of freedom near T c affect high px > 10 GeV jet flavor observables? 
Can an effective quasi-particle description be constructed that is consistent with lattice 
QCD thermodynamic data and simultaneously could reconcile the apparent inconsistency 
between the bulk “perfect fluid”, minimally viscous hydrodynamics and the success (mod¬ 
ulo V 2 ) of pQCD based jet quenching phenomenology? Can the combined set of soft plus 
hard observables be used to elucidate the mechanism of color confinement? The goal of 
this paper is to demonstrate in detail the existence of at least one possible model that 
involves partially suppressed cec together with emergent cmm effective degrees of freedom 
(the semi-Quark-Gluon-Monopole-Plasma model) that via the CUJET3.0 numerical frame¬ 
work allows us to constrain the q(E,T ) and p/s(T ) fields over a much wider range than 
previously possible. 

We begin by recalling a possible mechanism of color confinement based on a “dual 
superconductor” picture as proposed by Mandelstam, Nambu, Polyakov, ’t Hooft [39-42] 
It is understood that in type II superconductors the condensate of cooper pairs generates 
a “Meissner Effect” that repels magnetic field lines and squeezes monopole pair fields into 
flux tubes. In models possessing electric-magnetic duality, e.g. the Seiberg-Witten solution 
of the J\f = 2 supersymmetric gauge theory [43, 44], a “dual superconductor” phase does 
emerge in the strong coupling regime. Stable magnetic monopoles condensates can be 
generated leading to “Dual Meissner Effect” that forces the chromo-electric fields sourced 
by cec pairs to form flux tubes that gives rise to linearly rising potential and confines 
quark-antiquark pairs. 

Partially motivated by the “dual superconductivity” picture of color confinement, a 
magnetic scenario for the near T c QCD plasma was proposed in [45-50]. This scenario em¬ 
phasizes the change in chrorno degrees of freedom with the gauge coupling, and recasts the 
QCD phase diagram into electrically and magnetically dominated regimes. For example, 
focusing on increasing temperature (at zero baryonic density ), a particular temperature 
Te=m may be identified as a new phase boundary where the coupling strength of electric 


- 3 - 


(E) and magnetic (M) interactions are equal and satisfy the Dirac quantization condi¬ 
tion [51]. Interestingly, a novel post-confinement non-conformal non-perturbative region 
emerges at T c < T < Te=m ~ 1.4 T c where chromo magnetic monopoles (cmm) are the 
lightest degrees of freedom (DOFs) in the system while chromo electrically charged (cec) 
quasi-particles are strongly correlated and connected by flux tubes. Phenomenologically it 
has been qualitatively demonstrated that with the inclusion of such emergent monopoles 
near T c , the scatterings in both the soft and hard sectors are strongly enhanced and thus 
help explaining the small r\/s of the bulk sQGP as well as leading to significant V 2 of high 
Pt hadrons [45, 49, 52-55]. 

However, a quantitative and phenomenologically robust modeling framework for the 
QCD matter in the near T c regime has not been built previously. Such a framework has 
to couple the hard probes with the state-of-the-art bulk constrained viscous hydrodynamic 
(T(x, t), u^(x, t)) fields. It also needs to have a sophisticated implementation of the micro¬ 
scopic nonperturbative physics for both the cec and cmm degrees of freedom. For that, we 
need to constrain how the E and M quasi-particles DOF vary with temperature in a way 
that is consistent with available lattice QCD data [56-62] on thermodynamic pressure, 
entropy density, energy density as well as the Polyakov loop and quark susceptibilities. 
Specifically for the nonperturbative dynamics of the cec near T c , the “semi-QGP” model 
[63-68] was previously proposed to build in the color suppression effect in the deconfinement 
transition region. In that model the Polyakov loop is the relevant parameter that controls 
confinement/deconfinement of color electric charge DOFs. When temperature drops to¬ 
wards T c , the excitation of cec quarks and gluons are quenched by powers of the loop, 
resulting in a number of novel phenomenological effects. For example it was found in the 
semi-QGP that there is a mild enhancement and a strong suppression for the production 
of thermal dileptons and hard photons [69, 70]. 

The above considerations of the nonperturbative medium near T c have motivated us 
to propose and study a nonperturbative semi-Quark-Gluon-Monopole-Plasma (sQGMP) 
model for the “perfect chromo fluid” near the deconfinement transition range. In [71], we 
summarized the results of our CUJET3.0 extension of a pQCD based energy loss model 
CUJET2.0 discussed[38, 72] which integrates local parton energy loss over (2+l)D viscous 
hydrodynamic flows and and models jet medium interactions via the sQGMP quasi-particle 
model picture of the chromo structure of the fluid that include specific non-pertubative fea¬ 
tures related to confinement in the vicinity of T c . In CUJET3.0 all thermodynamic proper¬ 
ties are constrained by lattice QCD data. At very high temperature T>T C the model by 
construction would smoothly reduce to CUJET2.0 corresponding to a perturbative Hard 
Thermal Loop (HTL) picture of the QGP. As T approaches T c , the chromo-electric charge 
(cec) degrees of freedom are suppressed as powers of the Polyakov loop as in the semi-QGP 
framework while chromo-magnetic monopoles emerge to account for the total lattice QCD 
pressure or the total entropy density. Jet quenching observables of high px light hadrons’ 
and open heavy flavors’ Raa and at RHIC and LHC are computed using CUJET3.0 
and are shown to be in agreement with all data simultaneously withing present experi¬ 
mental uncertainties. This model therefore provides a semi-quantitative bridge between 
local equilibrium bulk “perfect fluidity” and high energy far from equilibrium jet transport 
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phenomena. 

Following up the condensed summary of the CUJET3.0 results reported in [71], we will 
present in this paper the theoretical details of the CUJET3.0 framework and concentrate 
on the robustness and consistency of its components as well on its phenomenological ap¬ 
plications that successfully account for a large set of leading-hadron-suppression data. We 
address a number of key questions concerning the theoretical robustness of the underlying 
sQGMP microscopic scenario and report new results that help in estimating theoretical 
uncertainties in our present understanding of jet energy loss phenomena via the sQGMP 
model. Such questions include: (1) Can an effective quasi-particle chromodynamic model 
be formulated with sufficient ab-initio lattice QCD and experimental data constraints to 
explain simultaneously both long and short wavelength observables in high energy A+A 
reactions? (2) How do theoretical uncertainties on Quark and Glue cec quasi-particle num¬ 
ber densities pq(T ) and Pg{T) and chrorno Magnetic monopole cmm density Pm(T) of the 
sQGMP propagate to the observables? (3) How do uncertainties in the chrorno electric 
and magnetic screening masses {pe,m{T ~ T c )) near T c effect the observables? (4) Can 
we constrain in the sQGMP model the temperature dependence of the jet quenching pa¬ 
rameter ( q(E,T )) well enough to predict the shear viscosity to entropy ratio ( r]/s(T )) by 
extrapolating down to E ~ 3 T thermal scales? (5) How does the effective path length 
dependence of light and heavy quark energy loss depend on the detailed cec and cmc com¬ 
position of sQGMP matter, and how do these compare to perturbative QCD HTL path 
length dependences? 

The organization of this paper is as follows: in section 2 , we briefly review the CU- 
JET2.0 pQCD energy loss kernel - the dynamical Djordjevic-Gyulassy-Levai-Vitev (DGLV) 
opacity series with multi-scale running strong couplings, the semi-QGP model, the magnetic 
scenario of the near T c QCD matter, as well as how the perturbative and non-perturbative 
ingredients are integrated in the CUJET3.0 framework. In section 3, we discuss two differ¬ 
ent scenarios for the deconfinement of color electric charge (cec) degrees of freedom near T c , 
and explore how the jet quenching observables computed from three schemes with different 
c m would vary in CUJET3.0 framework. We then investigate the influence of color com¬ 
position of the fluid on key transport properties, the q/T 3 and 77 /s, and propose schematic 
strategies for constraining 77 /s{T) from high px jet quenching data and vice versa. In 
section C, we present a systematic study of the path length dependence of parton energy 
loss in the sQGMP, including both the light quark and the heavy flavor. We summarize 
and propose possible future studies in section 4. We discuss in appendix A the effects of 
relativistic hydrodynamic flows fields on anisotropic jet suppressions in the sQGMP, and 
study in appendix B which single ingredient is most significant to generate the simultaneous 
satisfactory description of high px Raa and ^2 within CUJET3.0 framework. 

2 The CUJET3.0 framework 

The CUJET3.0 jet quenching framework generalizes the CUJET2.0 perturbative non- 
Abelian gluon bremsstrahlung kernel of the DGLV opacity expansion theory (section 2.1) 
of elastic radiative jet energy loss to incorporate lattice QCD constraints on the nonpertur- 
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bative semi-QGP suppression of color DOFs (section 2.2.1), and emergent nonperturbative 
chromo-magnetic monopole DOF near the crossover QCD temperature regime T ~ T c 
regime (section 2.2.2). CUJET3.0 incorporates lattice QCD data on the QCD pressure 
equation of state (EOS) P(T), nonperturbative chrorno electric and magnetic screening 
masses, m,E(T),mM(T ), ,and the Polyakov loop L(T) (section 2.3). In the following sub¬ 
sections, we present comprehensive discussions of the details of these ingredients and study 
robustness to systematic theoretical uncertainties 3.1 associated with “slow” quark liber¬ 
ation as suggested by L(T) data compared to “fast” quark liberation as suggested by the 
light quark susceptibility, y u (T), data. 


2.1 Perturbative QCD jet quenching: DGLV opacity expansion 


In the pQCD paradigm, radiative processes dominate the jet-medium interactions for a 
highly energetic parton. Based on different physical assumptions and approximations about 
the virtuality and branching of the hard parton, the nature of the medium through which 
the energetic parton propagates, and the kinematics for the interaction between medium 
and projectile parton [20], energy loss models like AMY [73-75], ASW [76-80], BDMPS- 
Z [81-83], Higher Twist [84-86], (D)GLV [17, 87-94] have been built and developed. In 
the dynamical DGLV opacity expansion theory [87, 89, 92, 95], the inclusive single gluon 
emission spectrum in the n = 1 opacity series with multi-scale running strong couplings 
[96] reads [38, 71, 72]: 


x E 


dN™= 1 

dXE 


18Cr 4 + Nj 
7 r 2 16 + 9IV/ 


x a s ( 


k 2 


dr p{ z)r(z) 

«2(ql) 


d 2 k ± 


2 ,, 2 / 


*+(!-*+) ■ 
— 2(k ± - q ± ) 


d A q- 


f E I 1 ' 


/i 2 (z) qi(qi + /!// 2 (z)) 
kj_ (kj_ - q±) 


1 — cos 


(k_L - q ±) 2 + X 2 (z) , 

2x4,E 


( k ± - 

q±) 2 


dx .| 

\x+) 

dXE 


( 2 . 1 ) 


Ce = 4/3 or 3 is the quadratic Casimir of the quark or gluon; the transverse coordinate 
of the hard parton is denoted by z = (xo + r cos <f>, yo + r sin </>; r); p{ z) and T(z) is the 
particle number density and the medium temperature in the local rest frame. In the 
presence of hydrodynamical 4-velocity fields Uj( z), boosting back to the lab frame, one 
should take into account a relativistic correction T(z) = [97, 98], where the flow 

4-velocity u^ = 7/(1, /3/) and null hard parton 4-velocity = (1 ,/?/)• E is the energy of 
the hard parton in the lab frame, kj_ (|kj_| < xeE ■ T(z)) and qj_ (|qj_| < 6T(z)E ■ T(z)) is 
the local transverse momentum of the radiated gluon and the local transverse momentum 
transfer respectively. The Debye screening mass p{z) is determined from solving the self- 
consistent equation 


p 2 (z) = sj 47ro: s (/t 2 (z))T(z)-y/l + Nf/6 (2.2) 

as in [99]; x 2 (z) = M 2 x 2 + +m 2 { z)(l—x + ) regulates the soft collinear divergences in the color 
antennae and controls the Landau-Pomeranchuk-Migdal (LPM) phase, the gluon plasmon 
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mass m g ( z) = /ep{ z)/\/2, and /e is the HTL chromo-electric deformation parameter, with 
/e = 1 by default [95]. The gluon fractional energy xe and fractional plus-momentum x + 
are connected by x + {xe) = xe [1 + y/l — (k±/xEE) 2 ]/2. 

In the CUJET2.0 model, Zakharov’s 1-loop pQCD running scheme is used [38, 100, 
101]. This running is cutoff in the infrared when the strong coupling strength reaches a 
maximum value a max for Q < Q m in- 


a s (Q 2 ) 


Oir 


47T 


k 91°g(Q 2 /Ag CD ) 


if Q < Qmin i 
if Q > Qmin • 


(2.3) 


where the minimum running scale Qmin is fixed by a max via Q m i n = Aq C d exp {2vr/9 a max } , 
with A qcd = 200 MeV. Note that the one power of a s originating from gluon radiation 
vertex runs with the virtuality fc^/[.x+(l — £+)], while the two powers of a s originating 
from parton-medium scatterings run with the transverse momentum exchange q\. 

For collisional processes, TG elastic energy loss formula [102-104] with Peigne-Peshier 
running coupling prescription [104] is used [38, 72]: 


= - C R n [a(n(z))a(QE(z)T(z)T(z))\ T{ z) 2 (l + ^ 


x log 


6T(z) V / F(z) 2 r(z) 2 - M 2 
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(2.4) 
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7r 


(4 + N f )T(i 


(2.5) 


respectively. Note that the calculation of the energy loss and the average number of colli¬ 
sions N c involves recursively solving the E( z) integral equation. 

In Eq. (2.1)(2.4)(2.5), the bulk evolution profiles (T(z), p(z),u^(z)) are generated from 
the VISH2+1 code [105-107] with MC-Glauber initial condition, to = 0.6 fm/c, s95p- 
PCE Equation of State (EOS), rj/s = 0.08, and Cooper-Frye freeze-out temperature 120 
MeV [108-113]. Event-averaged smooth profiles are embedded, and the path integrations 
J dr for jets initially produced at transverse coordinates (xo,</>) are cutoff at dynamical 
r(z(xo,</>,r))| Tmax = T cut = 160 MeV hypersurfaces [38]. 

Fluctuations about the mean radiative and elastic energy loss are taken into account in 
the following approximations: Poisson multiple gluon emission is assumed in the radiative 
sector; Gaussian fluctuations are assumed in the elastic sector. The total energy loss 
probability distribution is constructed from the convolution of the radiative and the elastic 
sector. In order to get the quenched leading hadron spectra at high px, this parton-level 
distribution is convoluted with LO pQCD pp spectra (for gluon and ligh quark) [114] 
or FONLL pp spectra (for charm and bottom) [115], folded over the MC-Glauber A+A 
initial hard scattering probability distribution [106, 116-119], and hadronized with parton 
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Figure 1. (Color online) (a) The measured invariant differential cross sections of inclusive charged 
particles in the mid-rapidity regime in p+p snn = 2.76 TeV collisions at the LHC (ALICE (red, 
[123]), ATLAS (blue, [124])), and those of neutral pions in p+p Snn = 200 GeV collisions at RHIC 
(PHENIX (orange, [125])), compared with CUJET pp references that are generated from LO pQCD 
calculations [114] with CTEQ5 Parton Distribution Functions (PDFs) and LO KKP Fragmentation 
Functions (FFs). (b) Data of prompt D° invariant differential cross sections (mid-rapidity) in p+p 
\JSnn = 2.76 TeV collisions at the LHC (ALICE (red, [126])) compared with CUJET pp references 
generated from FONLL calculations [115] with CTEQ6M PDFs and Peterson FFs. 


fragmentation functions (LO KKP [120] for gluon and light quark 3 , Peterson [122] for 
charm and bottom). 

The pp baselines embedded in CUJET are plotted in Fig. 1. For light hadrons, data 
from corresponding RHIC and LHC measurements [123-125] are well described by the LO 
pQCD results with a K-factor of 1.13 (fitting to the PHENIX datum at px = 6.24 GeV) 
at RHIC and of 1.76 (fitting to the ATLAS datum at pt = 6.26 GeV) at the LHC. In 
the heavy flavor sector, the FONLL results with a K-factor of 0.43 (fitting to the ALICE 
datum at pt = 7 GeV) are in agreements with data from ALICE measurements [126]. 
The most important message demonstrated by the comparison in Fig. 1, is that the px 
spectra shapes used in CUJET are consistent with RHIC and LHC data. The fitted LO 
K-factors (that aretheoretically predictable only with NLO and beyond) drop out of the 
single hadron observables Raa and V 2 which are only dependent on the spectra shapes and 
independent of absolute normalization. 

In the CUJET2.0 framework, only a = q, g (where a is the quasi-parton type, q stands 
for quark while g for gluon) HTL components are included. While in the CUJET3.0 

3 There are systematic studies that confront current NLO parton fragmentation functions with inclusive 
charged-particle spectra at hadron colliders and compare the varied consistencies between NLO FFs and 
data, cf. e.g. [121]. The light-hadron pp references in CUJET are generated from LO pQCD calculations 
with CTEQ5 PDFs plus LO KKP FFs. The results are consistent with available data, as shown in Fig. 1(a). 
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framework, a = m (m is the particle monopole) is also included and running coupling 
elastic and radiative energy loss are modified as we discuss in the next section to incorporate 
nonconformal nonperturbative lattice QCD data to constrain the composition and screening 
properties of the sQGMP. 

2.2 Nonperturbative QCD matter near T c 

In the vicinity of the QCD deconfinement transition temperature T c ~ A qcd, the strength 
of the chromo-electric coupling becomes very strong. Novel nonperturbative effects should 
enter and modify the properties of the QGP in this regime. Based on first-principle lattice 
QCD calculations and reasonable theoretical assumptions, in the CUJET3.0 framework, 
we model the near T c QCD matter as a semi-Quark-Gluon-Monopole-Plasma (sQGMP) 
that includes two ingredients with nonperturbative origins - the semi-QGP suppression of 
color electric DOFs and the emergent chromo magnetic monopoles. 

2.2.1 The semi-quark-gluon plasma model 

The semi-QGP model was developed and discussed in detail in [63-70]. It is constructed 
to describe QCD as temperature T — > T+, where both the naive perturbative methods 
and the hadronic models are not applicable. A main emphasis is on the “mismatch” 
seen from lattice data between the liberation of thermal excitations (as indicated by e.g. 
rapid increase of entropy density around T c ) and the liberation of “color” (as indicated by 
the rather slow increase of Polyakov loop toward fully deconfined limit). This observation 
indicates at a region above T c where significant nonperturbative suppression of color charge 
is still present. In the semi-QGP, how color is suppressed is quantified by the decrease of 
the expectation value of the Polyakov loop (£), which is the trace of a straight Wilson line 
in imaginary time. Properly normalized {£) is near unity in the perturbative QGP, but 
in the near T c regime, it is smaller than 1 (from lattice, {£) ~ 0.1 at T c ). This implies a 
non-trivial distribution for the eigenvalues of the Wilson line and a nontrivial background 
field for the time-like component of the gauge held Aq. In the presence of a nonzero Aq, 
as T —> T c + , the colored excitations are suppressed by powers of the Polyakov loop. 

Let us briefly review how this suppression works in the semi-QGP following [65, 66]. 
Adopting the double line basis for color factors, fundamental quarks carry a single index in 
the fundamental representation, a = 1, • • • ,N C , and adjoint gluons carry a pair of funda¬ 
mental indices, ab. In an SU(N C ) gauge theory, under mean held approximation, we take 
the temporal component of the gluon held to be a constant, diagonal, traceless matrix 

{A$) ab = 5 ab Q a /g (2.6) 

where g is the coupling constant for the SU(N c ) gauge theory. For the spatial components 
of the gluon vector potential, A t , there is no background held. In the Euclidean spacetime, 
the Wilson line in the temporal direction is 

( f 1/T 

L(x) = V exp ig j 


<1tAq{t , x) 


(2.7) 
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We neglect fluctuations in Aq to leading order in the coupling constant under mean field 
approximation. Gauge invariant Polyakov loops are traces of powers of the Wilson line: 


N c 

£n(Q) = (trL n )/N c = £ e in ^' T /N c . 

a= 1 


( 2 . 8 ) 


We define £ as the first Polyakov loop £ = £\. Physically, one can think of the Polyakov 
loop as measuring the extra free energy F which arises from adding a colored heavy quark 
to a thermal bath, {£) ~ exp(— F/T). In the perturbative QGP, all Q a, s vanish and £ n = 1. 
In the confined phase of a pure gauge theory, eigenvalues of Q are uniformly distributed on 
a circle of radius T, and £ n = 0 if n / [N C /2]N C . Dynamical quarks act as a background 
Z(N ) field, if they are present, there is no rigorous definition of a confined phase, and all 
Polyakov loops are nonzero at nonzero temperatures [65]. Lattice simulations find that £ 
is small {{£) ~ 0.1) in the phase transition regime for N c = 3 and Nf = 3 [66]. 

In the imaginary time formalism, the Euclidean four momentum is P ^ = (j>q, p), where 
Po is an even/odd multiple of nT for bosons/fermions. Expanding around the background 
field in Eq. (2.6), the 4-momentum of a quark becomes = (po + Q a , p), and the 4- 
momentum of a gluon becomes P® b = (po + Q ab , p) ( Q ab = Q a — Q b ). Each Q a is typically 
a non-integral multiple of 2irT, in the space of diagonal generators, it is like an imaginary 
chemical potential for color charges. To analytically continue from Euclidean to Minkowski 
spacetime, one continues the entire Euclidean energy to —iE, where E is a continuous 
energy variable. For quarks, the usual Fermi-Dirac statistical distribution function n(E ) 
becomes 


n, 


,.(E) = 


1 


e {E—iQ a )/T _|_ i 


EmP 1 ' 

3 = 1 


,—j{E—iQ a )/T 


(2.9) 


The first term represents the Boltzmann approximation to the quantum distribution func¬ 
tion, and is accompanied by exp (iQ a /T). Consider the trace of the quark propagator which 
enters e.g. in the calculation of the pressure at leading order. The sum is 


1 

N~c 




J](-l ) j+1 e ~ jE/T 

3=1 


1 

N~c 


trLi . 


( 2 . 10 ) 


Denote • - )q as the average over the Q distribution and an integration over the particles 
three momenta. At T —> T+ where £ is nonzero but small, 


(Y^na) Q ~N c T 3 £. ( 2 . 11 ) 


This means the distribution function for a quark field vanishes as a single power of the 
loop, i.e. n q ~ £. For gluons, the usual Bose-Einstein statistical distribution function n(E) 
becomes 


IT'ab(E') 


_ t _ = V p-j(E-i(Q a -Q b ))/T 

e (E—i(Q a —Q b ))/T l Z^ 

3=1 


( 2 . 12 ) 
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Consider summing over the color indices of the gluon propagator, to avoid taking the 
trace which is part of the 1/N C correction, one sums separately over a and b. Since 

e l< ^ a / T Ylb e ~ l ® b t T = trL trL 1, we have 

CtTn ab ) Q ~ N*T 3 f . (2.13) 

ab 

This means in the near T c regime, the density of gluons vanishes as the square of the loop, 
i.e. n g ~ £ 2 . In the perturbative regime, the density of massless fields is necessarily a pure 
number times T' 3 therefore such a suppression is not present. Note that we always perform 
a global color rotation to enforce that the expectation value of the Polyakov loop l is real. 

2.2.2 The magnetic scenario for the near T c QCD plasma 

The magnetic scenario for the near T c QCD matter was proposed and discussed in [45-50], 
and there have since been extensive studies of the magnetic component of the plasma using 
different approaches [52, 127-129]. In this scenario, the QGP not too far above T c contains 
not only electrically charged quasi-particles (cec), quarks and gluons, but also magnetically 
charged quasiparticles (cmc), monopoles and dyons. 

This approach is different from many traditional discussions, which focus on the ther¬ 
modynamic transition and divide the temperature regimes into the hadronic phase at 
T < T c and the QGP phase at T > T c . Rather, the emphasis is on the competition 
between EQPs and MQPs, based on which one may divide the phases of QCD matter into 
the magnetically dominated region at T < Te=m ~ 1.4T C and the electrically dominated 
region at T > Te=m • This picture is largely motivated by analogy with electric-magnetic 
duality in supersymmetric Yang-Mills theories. The key aspect of the physics involved is 
the coupling strength of the electric (e) and magnetic (g) interaction, which can lead to 
different dominance of dynamical degrees of freedom in different regimes. Under Dirac 
quantization condition [51], e ■ g = n/2, and magnetic objects are in the adjoint color rep¬ 
resentation if n = 2. In a so-called E/M-equilibrium region, the couplings are equal, i.e. 
e = g, densities as well as masses of both EQPs and MQPs are comparable. Then depend¬ 
ing on the change of these couplings in different physical regimes, the “balance” between 
E and M sectors would shift one way or the other, giving rise to distinctive phases. 

Let us start with the QGP at very high temperature T T c where the electric 
coupling is weak. This regime is well described by perturbative EQPs with small quark and 
gluon effective masses. The monopoles in this case are heavy, dilute and strongly coupled, 
but they play a minor role and contribute little to the overall bulk properties. They do 
manifest themselves through nonperturbative contributions to certain observables at the 
soft magnetic scale. 

On the other hand, as T goes down and approaches the confinement transition T —> T c , 
the converse is expected to happen: the electric coupling becomes very strong and the 
EQPs, i.e. quarks and gluons, are getting heavier and gradually suppressed dynami¬ 
cally. The emergent MQPs gradually become light, abundant, and dominate the system 
at T < Te=m ■ With further decrease of temperature toward T c the thermal monopoles 
will eventually reach the Bose-Einstein condensation, forming a ’t Hooft-Mandelstam “dual 
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superconductor” [39, 42] that enforces color confinement. In the post-confinement region 
at T c < T < Te=m EQPs are still strongly correlated and connected by the electric flux 
tubes, but MQPs are the dominant DOFs and they serve as an effective description of the 
strong nonperturbative gauge dynamics. In [48], the authors showed that gauge theory 
monopoles in a deconfined phase behave as magnetic charges in a Coulomb plasmas. At 
T ~ 1.3T C where lattice potentials indicate flux tubes dissolve, an estimate of total den¬ 
sity of magnetic quasi-particles is n cmc ~ 4.4 — 6.6 fm 3 [46]. From an analysis of the lattice 
monopole-(anti)monopole correlators, they showed that the temperature dependence of the 
magnetic couplings in gauge theories is indeed the inverse of the electric one as per the 
electric-magnetic duality arguments. More specifically, the magnetic part of the QGP at 
T ~ 1 — 3T C possesses an effective plasma parameter in the “good liquid” domain, thus in 
consistency with the “nearly perfect liquid” property observed at RHIC and LHC. 

2.3 Jet suppression in semi-Quark-Gluon-Monopole-Plasmas 

Having discussed the foundations of the sQGMP, let us integrate it into the jet energy 
loss kernel in section 2.1. The critical component in Eq. (2.1) is the 1-HTL dynamical 
scattering potential, 



(2.14) 


Since the sQGMP contains both chromo electrically charged quasi-particles (cec) and 
chrorno magnetically charged quasiparticles (cmc), when jets propagate through the medium 
near T c , scattering channels of E + E and E + M exist simultaneously. One way to gen¬ 
eralize Eq. (2.14) is to symmetrize it with respect to the E and M components of the 
kernel based on demanding electric-magnetic duality as illustrated in e.g. the celebrated 
Seiberg-Witten solution of the AT = 2 super-Yang-Mills theory. This leads to the following 
modified form of the kernel: 



dN 


(2.15) 


x —— oc .. 
dx 


Where a s = a®, and cue • olm = 1 at any scale by Dirac quantization condition [48, 51]. 
The total quasi-particle number density p is divided into EQPs with fraction \t = Pe/P 
and MQPs with fraction 1 — \T = Pm/P- The parameter and /m is defined via /e = 
Pe/P and /m = Pm/P, with pe and pm being the electric and magnetic screening mass 
respectively. We emphasize that Eq.(2.15) is a nonperturbative sQGMP model ansatz that 
differs substantially from other generalization of HTL, see e.g. [95, 130]. 

To determine XTi one notices that: (1) when temperature is high, xt should reach 
unity, i.e. xt(T 3> T c ) —> 1; (2) in the vicinity of the regime T ~ (1—3)T C , the renormalized 
expectation value of the Polyakov loop L (let us redefine L = t = (trV exp{ ig f ( ] ' T dr A)}} / N c ) 
deviates significantly from unity, implying the suppression ~ L for quarks and L 2 for 

gluons in the semi-QGP model [64, 66, 68, 131]. In the regime the quark and gluon density 
drop much faster than the thermodynamic quantities. This points to “missing” DOFs, 
in the magnetic scenario [45, 48], they are identified as chromo-magnetic monopoles who 
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emerge in gauge theories at strong coupling and are thermal excitations of the vacuum mag¬ 
netic condensate as in the “dual superconductivity” picture of color confinement [132-134]. 
For the cec component fraction, we use the semi-QGP ansatz: 

Xt{T) = c q L(T) + c g L 2 (T) . (2.16) 


For the respective fraction of quarks and gluons, where we take the Stefan-Boltzmann 
(SB) fraction coefficients, c q = (10.5Nf)/(10.5Nf + 16) and c g = 16/(10. 5Nf + 16). To be 
consistent with lattice data, we fit the temperature dependent Polyakov loop L(T) (T in 
GeV) with 


L(T) 


1 + *Tanh[7.69(T - 0.0726)] 


(2.17) 


Eq. (2.17) adequately fits both the HotQCD [59] and Wuppertal-Budapest [60] Collab¬ 
oration results, c.f. Fig. 2(a). With xt and (1 — xt), Pe/P and pm/P are completely 
fixed. 

To specify the electric and magnetic screening mass (pe,m = Je,m p), we recall that 
at very high temperature, one expects (1) f E 1, i.e. pe ~ gT from HTL results and (2) 
/m ~ 55 i-e. pm ~ <? 2 T, from magnetic scaling in dimensional reduction. Assuming E-M 
duality, the screening masses are expected to scale as 


Pe,m 


~ che,m Pe,m/ t ■ 


(2.18) 


The extrapolation to lower temperature thus gives 

Pe ~ ole{xtp)/T ~ xtP 2 , ( 2 -19) 

and we expect the electric screening mass to be suppressed as y/xriT) in the near T c 
regime but approach the HTL p(T) at high T limit. Regarding the magnetic screening 
mass, since we have um ~ (c) 3 following the magnetic scaling, then 

Pm ~ « M ( a E Tf/T ~ a 2 E T 2 ~ g 2 g 2 T 2 ~ g 2 p 2 . (2.20) 


This prescription is supported by lattice data [58]. Therefore, we assume the following 
local temperature dependent screening masses in the CUJET3.0 model: 

f E (T( z)) = Vxt(T(z)) , f M (T(z)) = c m g(T(z)). (2.21) 


To be consistent with previous treatments in Eq. (2.1) and (2.2), the local electric “cou¬ 
pling” g(T(z)) is defined via 

s(TM) = Am?™)) = r7=f (2.22) 

T(z) v /l + Af / /6 

Note that c m is a constant parameter that can be constrained by lattice data on the 
magnetic screening. Fig. 2(b) illustrates the agreement between this prescription of pe,m 
and lattice extracted values [58]. 
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Finally, in the CUJET3.0 framework, the energy loss kernel Eq. (2.1) is generalized to 


x E 


dN™= 1 

dXE 


x 

X 


X 


18Cjj 4 + Nf 
tt 2 16 + 9 Nf 


d\ 


a 


l(qi) (/. 


dr p(z)r(z) J d 2 k±a s ( 
fl ft, ft ft) 


X+0--X+) 

2 f 2 .,2 


) 


2 + £i£ijp) ^ + + (1 _ XT ) 


(ql + /|M 2 (z))( qi + /> 2 (Z)) 

k_L (kj_ - qjJ 


-2(kj_ - qj_) _ 

(kj_ - q_i) 2 + X 2 (z) |_k 2 + x 2 (z) (kj_ - q±) 2 + x 2 (z) 
1 I (kj. -q±) 2 + X 2 (z)_ 

1 ^ cos| - 2EX - T 


( k _L - 

q±) 2 

1 ( XE \ 

dx .| 

\*+/ 

dXE 


(2.23) 


where xt and /e,m follows Eq. (2.16) and Eq. (2.21). We note that in the temperature 
range T ~ T c , the coupling a s becomes non-perturbative [45, 48, 100, 135]. Analysis of 
lattice data [48] suggests the following thermal running coupling form: 


a s (Q 2 ) = - 

1 + 


OLq 



(2.24) 


with T c = 160 MeV. Note that at large Q 2 , Eq. (2.24) converges to vacuum running 
a s(Q 2 ) = 9 iog(QVA 2 ) ; while at Q = a s{T 2 ) = a c . 


3 Liberation schemes for color degrees of freedom 

As discussed in section 2.2, within the semi-QGP model, the expectation value of the 
Polyakov loop L (note that we have redefined L = l in section 2.3, and we use this notation 
for the rest of this paper) is the only relevant parameter for the confinement/deconfinement 
transition, upon proper renormalizations, L serves as a suppression factor for the colored 
excitations as T —> T+. However, it is questionable whether or not L is an order parameter 
for the phase transition. Besides the fact that lattice calculations point to a L ~ 0.1 at T c , 
Eq. (2.6) (2.7)(2.8) also indicate that to a certain degree the loop physically characterizes 
the free energy of an infinitely massive static quark. Since (1) in the perturbative QGP 
phase dynamical light quarks dominate the medium transport properties; and (2) to boost 
V 2 in line with data, a strongly enhanced jet scattering near T c makes decisive contributions 
[71]; then, the nonperturbative property of the sQGMP, in particular, the rate at which 
fractional chromo-electric DOFs are liberated (defined as rd(T ) = dxr/dT) in the near 
T c regime will play a significant role in computing jet quenching observables within the 
CUJET3.0 framework and should be studied more systematically. 

3.1 Polyakov loop versus quark number susceptibility 

Another useful measure of the nonperturbative suppression of the color electric DOF is 
provided by the quark number susceptibilities [136-139]. Such susceptibilities quantify 
the quark number fluctuations that can be obtained from the QCD partition function at 
vanishing chemical potentials. Denote u, d, s as up, down, strange quark whose numbers 
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are conserved charges in QCD. Starting from the pressure, 



(3.1) 


moments of charge fluctuations are defined as follows, 



Qi+j+kpjrpA 


(3.2) 


difiu/Tydina/Tydfjis/T)* • 


Concentrate on the quadratic fluctuation, 



(3.3) 


And X 2 ’ d ’ S i s the diagonal susceptibility of u, d, s quark number density. Singlet suscep¬ 


tibilities of other conserved charges in QCD such as baryon number B, strangeness S and 


electric charge Q can be obtained from the above quark number susceptibilities [62], 

The diagonal susceptibility is proposed as part of the order parameter for chiral sym¬ 


metry breaking/restoration in [136]. Considering a gas of free quarks, if the quark mass 
m is small, then X 2 is expected to be large since it is relatively easy to create an addi¬ 
tional quark. For instances, if m <C T, then in the continuum limit, X 2 ~ NfT. If m is 
large, then it will be difficult to create a quark or antiquark, the susceptibilities will be 
suppressed by exp(— m/T). Realistically, in the high T phase, though strongly interacting, 
if the fundamental excitations of the system are low-m objects with the quantum numbers 
of quarks, then X 2 is still expected to be large. Meanwhile, in the low T phase, X 2 will be 
small since quarks are confined and the nonzero quark number states have large masses. 
Thus in the chirally symmetric phase, the quark number susceptibility is large, which is 
consistent with a plasma of light quarks; while in the chiral symmetry broken phase, the 
quark number susceptibility is small, as expected from quark confinement. It however may 
be noted that in the parton-hadron boundary regime, various bound states like baryons and 
mesons (and even other exotic composite objects) carry conserved charges and contribute 
to the susceptibilities. As previous studies have shown [140-142], such contributions are 
mostly important for the higher-order susceptibilities as well as for the off-diagonal ones. 
The leading order diagonal susceptibilities could therefore serve as a reasonable measure 
for the counting of quark degrees of freedom in the plasma. 

Therefore, besides interpolating the renormalized lattice Polyakov loop as in Eq. (2.17), 
we parametrize the lattice diagonal susceptibility of u quark number density as 


X%(T) = 0.91 x i {1 + Tanh[15.65(T - 0.0607)]} 


(3.4) 


Where T is the temperature in the unit of GeV. Note that at extremely high temperature, 
the X 2 X) is not unity, so we renormalize the susceptibility by its value at T —> oo and 
define a new quantity XaX) as 


X2X) = \i 1 + Tanh[15.65(T - 0.0607)]} 


(3.5) 
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Figure 2. (Color online) (a) Renormalized Polyakov loop L(T ) (blue circle: [59], green square: 
[60]) and diagonal susceptibility of light quark number density x^T) (red diamond: [62]) computed 
from lattice QCD, fitted with the parametrization of Eq. (2.17) and (3.4). The inset shows the 
the density fraction of color electric DOFs (red, xt = Pe/p) and color magnetic DOFs (blue, 
1 — Xt = Pm/p) within the liberation scheme Xt (solid) and Xt (dashed), in the temperature range 
T ~ 0.6 — 6.0 T c , where T c = 160 MeV. Notice that in Xti Pe ~ Pm at T ~ 1.7 T c ; in Xti Pe ~ Pm 
at T ~ 1.1T C ; and these temperatures are where rd(T) = dxr/dT should peak in Xt’ U - 0 3 ) The 
dimensionless electric (red) and magnetic (blue) screening mass Pe,m/T in the CUJET3.0 model 
i.e. Eq. (2.21), for scheme (i) (3.7), (ii) (3.8), and (iii) (3.9), compared with results from Lattice 
QCD [58]. Note that the a c and c m parameters in (i)(ii) (iii) are chosen such that the high-pr 
reference Raa datum can be well-fitted, c.f. Fig. 3(a). Note that for (i) and (ii), pe/T ~ Pm/T 
at around the same temperature, i.e. T ~ 1.5 — 1.6T C , while (iii)’s pe/T and pm/T intersect at 
T ~ 1.1 T c . In the near T c regime, (i) and (ii)’s pe — Pm are approximately identical, both are less 
than (iii)’s. 


The ^(T) Plays a similar role as properly renormalized L for quark DOFs. Let us denote 
the original liberation scheme, c.f. Eq. (2.16), that follows the power law of the Polyakov 
loop as in the semi-QGP model, as Xt (Xt = XT in Eq. (2.16)); and the new deconfinement 
scheme where the diagonal susceptibility of light quark number density dominates the 
transition, as Xt (Xt = Pe/p) : 

Xt = c qX2+ c g 1/2 ■ ( 3 - 6 ) 

Note that in this scheme, the magnetically charged quasi-particles, i.e. chromo-magnetic 
monopoles, consist a density fraction of 1 — Xt = Pm/ P- 

Fig. 2(a) shows the lattice QCD data on renormalized Polyakov loop and diagonal light 
quark susceptibility from the HotQCD [59] and Wuppertal-Budapest [60, 62] Collaboration, 
as well as the parametrization fit Eq. (2.17) and Eq. (3.4). The E and M quasi-particle 
density fraction in both the Xt anc ^ Xt scheme are plotted in the inset of Fig. 2(a). Note 
that pe/P = Xt* an d Pm/P = 1 — Xt’ U - The two different schemes, for the rate of “quark 
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liberation”, with \t the “slow” and Xt the “fast”, provide useful estimates of theoretical 
systematic uncertainties associated with the quark component of the sQGMP model. Note 
that in the inset of Fig. 2(a) the pe and pm are equal at T ~ 1.1 T c for Xt scheme while at 
T ~ 1.7 T c for Xt scheme: these temperatures are where r^(T) = / dT should peak. 

As suggested in Eq. (2.21), any change in xt alters the electric screening mass pe, 
and one expects shifts in the magnetic screening correspondingly under electric-magnetic 
duality. In Fig. 2(b), lattice data of the electric and magnetic screening mass are compared 
with the CUJET3.0 results in three schemes: 


(i) 

a c = 0.95, c m = 0.3, Xt ! 

(3.7) 

(ii) 

a c = 0.95 ,c m = 0.4, Xt! 

(3.8) 

(hi) 

a c = 1.33, Cm = 0.3, Xt ■ 

(3.9) 


Note that the ( a c , c m ) parameters are chosen such that the single reference datum R h AA (pT = 
12.5GeV) ~ 0.3 at LHC is well-fitted, c.f. Fig. 3(a). Implicitly, Xt and Xt I s determined 
by Eq. (2.16) and Eq. (3.6) respectively. All three schemes are in reasonable agreements 
with the lattice data. However, to be more careful, (i)(ii)(iii)’s pe,m do behave differently 
as temperature varies. (i)’s and (ii)’s pe/T and pm/T intersect at approximately the 
same temperature, i.e. T ~ 1.5 — 1.6T C . Meanwhile, (iii)’s pe/T and pm/T intersect at 
T ~ 1.1 T c , which temperature overlaps approximately with the T where Pe(T) = Pm(T ) 
in Xt ■ Furthermore, in the near T c regime, (i)’s and (ii)’s pe — Pm are approximately 
equal, and both are less than (iii)’s. 

Let us move to the predictions of jet quenching observables in scheme (i)(ii)(iii) within 
the CUJET3.0 framework, in particular, leading light hadrons’ and open heavy flavors’ 
Raa(pt > 8GeV) (c.f. Eq. (1.1)) and v 2 (p T > 8GeV) (c.f. Eq. (1.2)) at RHIC and LHC 
semi-peripheral A-(-A collisions are of interests. The results and corresponding data are 
plotted in Fig. 3. 

For high pr light hadrons (e° at RHIC, at LHC), Fig. 3(a) shows that all three 
schemes can describe the Raa(pt) data at RHIC (8 < pt < 20 GeV) and LHC (8 < pt < 
100 GeV) simultaneously, but only (i) and (ii) are compatible with the high pt single 7r° and 
/i ± ’ s v 2 as illustrated in Fig. 3(c). Since (i) and (ii) have different liberation schemes hence 
s, different absolute values of pe,m , but approximately the same pe — Pm near T Cl this 
observation indicates that in boosting the 7r° or h^'s azimuthal elliptical harmonics to be in 
line with data, the difference between pe and pu in the near T c regime plays a critical role. 
Notice that as suggested in the magnetic scenario, when cooling down to pass T ~ 1 — 2T c , 
the lightest hence the dominant DOFs in the medium shift from EQPs to MQPs, and 
the color screening mass is one of the indicators of this transition [45]. The fact that 
(i)(ii)’s Pe(T) ~ Pm{T) generates a larger v 2 than (iii)’s Pe(T) — Pm(T) who has a larger 
value and a lower zero point temperature To (To is defined as pe{Tq) = pm(Tq)) implies 
that, v 2 is sensitive to how the relative value of pe and pm inverses near T c - the higher 
the To, the longer the path length that jets interact with the strongly coupled monopole 
dominated medium at later time of the evolution, the larger the azimuthal anisotropy. A 
further comment is on the absolute values of pe,m hi (i) (ii), clear (ii) ’s are larger. This 
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Figure 3. (Color online) (a) Neutral pion (7 r°, brown) and charge particle (ft , blue)’s 
Raa(pt > 8GeV) in Au+Au 200GeV and Pb+Pb 2.76TeV 20-30% collisions, computed from CU- 
JET3.0 with scheme (i) (3.7) a c =0.95, c m =0.3, Xt (solid), (ii) (3.8) a c =0.95, c m =0.4, Xt (dashed) 
and (iii) (3.9) a c =1.33, c m =0.3, Xt (dotdashed), compared with corresponding RHIC [1-4] and 
LHC [5-9] measurements. The 7r° and h ±, s V 2 (pr > 8GeV) are plotted in (c). (b) The CUJET3.0 
results of D meson (red) and B meson (green)’s Raa(pt > 8GeV) at LHC in (i)(ii)(iii) compared 
with available data [10-12], The D and B’s V 2 {pt > 8GeV) are plotted in (d). Note that the 
parameters ( a c ,c m ) are determined in scheme (i)(ii)(iii) such that the reference datum at LHC 
R-aa(pt = 12.5GeV) « 0.3 as well as the lattice pm ([58], c.f. Fig. 2(b)) are reasonably fitted. 
Results from all three schemes are compatible with light hadron (LH)’s Raa] while for LH’s V 2 , 
(i) and (ii) can generate reasonable agreements with data, but (iii) underestimates the V 2 - For 
open heavy flavors (HF), (ii) and (iii) have similar Raa predictions, both differ from (i); while 
for HF’s i> 2 , (i), (ii), and (iii)’s prediction are all different. Such differences in the predictions for 
jet quenching observables from (i)(ii) (iii) suggest that data on high pt Raa and V 2 can impose 
stringent constraints on the nonperturbative properties of the medium near T c . 


is necessary because after T drops lower than To, MQPs dominates, (i)’s pm is denser 
than (ii)’s, to get to the same magnitude of overall leading hadron suppressions, (ii) should 
possess larger color screening masses, c.f. Eq. (2.15). 

For open heavy flavors, specifically, high pp D and B mesons, Fig. 3(b) shows their 
respective Raa at LHC Pb+Pb ^snn = 2.76 TeV collisions, computed from CUJET3.0 
with scheme (i) (ii) (iii). Differ from the light hadrons’ Raa where (i) (ii) (iii) have almost 
identical predictions, for D and B mesons, (ii) and (iii)’s Raa overlap, both of which distin¬ 
guish from (i)’s. One notices that (ii) and (iii) have different pe,m , but the same liberation 
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scheme Xx- w hich is distinct from (i)’s Xx- This implies the open heavy flavor’s high 
Pt Raa is critically influenced by the rate at which chromo-electric DOFs are deconfined 
{rd = dxr/dT). This connection is intrinsically embedded in the CUJET3.0 framework 
because the heavier masses induced dead cone effects shuffled the weights of the scattering 
potential (2.15) and the rest of the Eq. (2.23) in such a way that the total induced radiation 
became more sensitive to the deconfinement scheme rather than screening masses He M 
for jet quenching in sQGMP. On the other hand, Fig. 3(d) shows the prediction of open 
heavy flavor’s V 2 (px > 8GeV) at LHC in CUJET3.0. The U 2 ’s are all different in scheme 
(i), (ii), and (iii). This suggests the open charm and beauty’s V 2 {pt) act as good probes 
of the nonperturbative (rd, He, Hm) near T c . Let us mention in passing that for the heavy 
quark dynamics in low px region, the sQGMP model also expects a strong temperature 
dependence of their in-medium diffusion coefficients (specifically a “dip” near T c ), which 
has recently been shown as an essential ingredient toward a simultaneous description of 
their Raa and V 2 [143, 144], 

If one views the above physical connections from a different angle, the set of RHIC and 
LHC heavy-ion collision data on high px light hadron (LH) and open heavy flavor (HF)’s 
Raa and V 2 will provide stringent limits on the nonperturbative properties of the QCD 
matter near T c in the CUJET3.0 framework. Specifically, after parameters in the model 
have been fixed by light quark’s Raa data, the rate at which color DOFs are deconfined 
rd and the color screening masses he,m can be stringently constrained: (1) light quark’s 
V 2 regulates he(T) — hm(T) near T c ; (2) heavy quark’s Raa determines r^(T); (3) HF’s 
V 2 distinguishes all rd(T), he{T ) and hm(T). 


3.2 Jet transport coefficient and shear viscosity 

As discussed above, the high px Raa and V 2 data of light and heavy quarks can provide 
stringent constraints on values of the nonperturbative (rd, He, Hm) near T c . It is of great 
interests to further compare how the jet and bulk transport properties differ in varied 
schemes (i) (3.7), (ii) (3.8), and (iii) (3.9). This will pave the way for clarifying the 
temperature dependence of jet quenching and shear viscous transport properties based on 
available high px data in high-energy A+A collisions. 

The jet transport coefficient q characterizes the averaged transverse momentum trans¬ 
fer squared per mean free path [22], Here let us first calculate the q for a quark jet (in 
the fundamental representation F) with initial energy E, in the same way as the previous 
CUJET3.0 computation in [71], via the following: 


rGET 


q F (E,T) = 


dqi- 


2vrq2 


-AT) 


J 0 (qi+/i/i 2 )(qi+/ A V)' 

X [Aqqfq + Cqgfg) a s ) + C qm (l ~ fq~ fg)] ■ 

where p(T ) is the total number density, connected to the lattice pressure p(T) via 


(3.10) 


p(T)=ip(T)/T , (3.11) 

with £ = [90C(3)(16 + 9A(y)]/[7r 4 (16 + 10.5Aj)] = 1.012 for an N c = 3, Nf = 2.5 ideal gas. 
The parameters f q , g are quasi-parton density fractions of quark (q) or gluon (g) type, in 
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the Xt an d Xt scheme, they are respectively 

f q = c q L(T), fg = c g L(T)\ if XT ; 

f q = c g x!(T), fg = CgL(T) 2 , if XT ■ (3.12) 


The c qt g and L(T ) are the same as in Eq. (2.16) and (2.17). The magnetically charged 
quasi-particle density fraction is hence f m (T) = 1 — f q (T) — f g {T). The color factors in 
Eq.(3.10) are given by 


C qq 

C q9 


4 

g’ l — Cmm — 

Cgq — Cgm — mg 


9 

4’ 

= 1 . 


(3.13) 


In the CUJET3.0 framework, once the jet transport coefficient q has been computed, 
the shear viscosity to entropy density ratio ij/s can be calculated based on kinetic theory 
in a weakly-coupled quasi-particle picture, as proposed in [30, 32, 146]. An estimate of p/s 
can be derived as 


v/s 


1 4 
s 15 


Pa{p)a^ a 


('ST' [ {Sab)/2 J 2 4 q]_ da ab 

5 , ^Pa/q a (T, E = 3T) . 


4 T 
5s 

18T 3 


-l 


(3.14) 


Note that the q(T,E ) is extrapolated down to thermal energy scales E ~ 3 T. The p a (T) 
is the quasi-parton density of type a = q,g,m. The mean thermal Mandelstam variable 
(Sab) 18 T 2 . The entropy density s(T) is interpolated from lattice calculations [59]. The 
q a=q = qp is calculated as in Eq. (3.10). The q a =g,m for are computed via 


q a = g (E, T) 

q a = m (E, T) 


[ btjl 27rq 2 

Jo dq± (q± + f 2 E p 2 )( qi + fhp 2 ) P ^ 

X [(C gq f q + Cggfg)a 2 s (ci 2 ± ) +C gm (l- fg- fg)\ , 
j-OET 27rq 2 

Jo dq± (ql + P e p 2 )( qi + fhp 2 ) P ^ 

X \(C’m q f q + C mg fg)+C mm (l — f q —fg)/a 2 ((P 1 f)] . 


(3.15) 


(3.16) 


Clearly the q/s of the system is dominated by the ingredient which has the largest p a /q a - 
Fig. 4(a)(b)(c) shows the dimensionless jet transport coefficient q/T 3 in CUJET3.0 for 
a quark jet with initial energy E = 20, 10, 2 GeV respectively, compared with the CUJET2.0 
result [38] and the AdS/CFT limit [145]. The q in scheme (i)(ii)(iii) are plotted. Compared 
with (i) which has xJf a sfo w ” quark liberation, the q in (ii) and (iii) which have Xt a fast” 
quark liberation scheme are significantly smaller. This may be understood as follows: in the 
Xt scheme, as temperature rises, the chromo-electric DOFs are excited faster than those 
in the Xt scheme, and leads to a smaller fraction of magnetically charged quasi-particles 
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Figure 4. (Color online) The temperature dependence of the dimensionless jet transport coefficient 
q/T 3 for a light quark jet (F) with initial energy E = (a) 20GeV, (b) lOGeV, (c) 2GeV in the 
CUJET3.0 framework (Red) with the three schemes: (i) (3.7) (solid), (ii) (3.8) (dashed), and (iii) 
(3.9) (dotdashed). The CUJET2.0 </f/T 3 with (a max ,fE, $m) = (0.39,1, 0) (Blue) and the J\f = 4 
Super Yang-Mills (SYM) qsym/T 3 = [145] with’t Hooft coupling A = 47 t (Black) are 

plotted for comparisons. The insets show the absolute in CUJET3.0, CUJET2.0 and SYM. Note 
that the (ii) and (iii) ’s q are similar, and both are smaller than (i)’s. (d) The shear viscosity to 
entropy density ratio 77 /s estimated in the kinetic theory using the q extrapolation Eq. (3.14) in 
CUJET3.0 with scheme (i) (solid) (ii) (dashed) (iii) (dotdashed), for quasi-patron type q (quark, 
red), g (gluon, green), and m (monopole, blue). The total r\/s is plotted with black curves. The 
inset shows quasi-particle number density fraction of q, g, m in the liberation scheme \t an d Xt- 
Note that in the near T c regime, in the \t scheme, the total i]/s is dominated by q, while in the 
Xf “slow” quark liberation scheme the total 77 /s is dominated by m. In addition, there is a clear 
77 /s ~ 0.12 minimum at T ~ 210 MeV in (ii) and (iii) which utilize the same \t “fast” quark 
liberation scheme. This ( rj/s) m i n is larger and phenomenologically more favorable than that in the 
“slow” quark liberation scheme. 


in the near T c regime. Since chromo-magnetic monopoles are the key contributors to the 
enhancement of jet opacity, c.f. appendix B, less monopoles thus result in a diminished q 
in the Xt scheme compared with the one in the Xt scheme. Interestingly, (ii) and (iii)’s q 
still get close to the Super Yang-Mills (SYM) limit near T c . Note that (ii) and (iii)’s he,m 
behave very differently as shown in Fig. 2(b), then a crucial observation one can draw is 
that among the nonperturbative ( rH e,Hm), the high energy jet transport property which 
is determined by the quenching parameter q(T), is sensitive to r^, he. the rate at which 
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confined colors are excited near T c . Apart from such sensitivity, the near-T c enhancement 
of jet-medium interaction is a very robust feature in all schemes and is deeply tied with 
the nonconformal, nonperturbative dynamics near the transition temperature. It may be 
worth mentioning that a number of studies based on holographic QCD models [147, 148] 
that build in near-T c conformal-breaking effects, have universally found similar near-T c 
enhancement of the q/T 3 . 

A surprising aspect of the comparison between CUJET2.0 and CUJET3.0 in Fig. 4, 
is that both models describe equally well the azimuthally averaged Raa data (which char¬ 
acterize overall suppression), in spite of their rather different functional forms of qp(T). 
The two models differ from each other most significantly in the near-Tc regime: the q-p/T :i 
of CUJET3.0 is much larger than that of CUJET2.0 for T ~ (1 — 2 )T C . Above T ~ 3T C 
and by T ~ 6 T C , the 3 of CUEJT2.0 is ~ 50% larger than that of CUJET3.0. The 
overall energy loss is controlled by qF (rather than qp/T s ) and therefore more dominated 
by contributions from the high temperature QGP in the fireball. This explains why both 
CUJET2.0 and CUJET3.0 are able to fit the overall Raa- The azimuthal anisotropy V 2 , on 
the other hand, is more sensitive to the late time contributions to energy loss coming from 
the lower temperature T ~ (1 — 2)T C part of the fireball. While CUJET2.0 fails to describe 
V 2 , the CUJET3.0 successfully describes the V 2 data precisely by virtue of the strongly 
enhanced near-T c contributions due to the emergent color magnetic monopole degrees of 
freedom in the sQGMP. The contrast between the CUJET2.0 and CUJET3.0 demonstrates 
again the importance of simultaneous descriptions for both Raa and V 2 data in order to 
differentiate energy loss models. 

The fact that remarkably different q(T ) dependence could be consistent with 
the same Raa data, demonstrates clearly the inadequacy of focusing on the jet 
path averaged quantity {q) as the only relevant medium property to characterize 
jet energy loss. Evidently while the {q) captures the important transverse “kick” factor, 
there are other essential factors like the actual chromo electric and magnetic composition 
of the plasma, the screening masses and the running couplings at multiple scales which all 
strongly influence jet energy loss and imprint their effects beyond just in the {q}. 

In Fig. 4(d) the q extrapolated r\/s in scheme (i)(ii)(iii) following Eq. (3.14) are plotted. 
Note that when T < T c , as T keeps cooling down r\js rises due to the hadron resonance 
gas (HRG) contributions, as computed in [149-152], In terms of the total q/s near T c , 
both (ii) and (iii) has a clear minimum ( r//s) m i n ~ 0.12 at T ~ 210 MeV, while (i) has 
a ( 77 / s)mi n ~ 0.02 which is under the 77 /s ~ 0.08 quantum bound. This suggests (1) the 
liberation scheme of color DOFs, i.e. r^, significantly influences the lower bound of 77 /s; 
(2) ( r\/s)min is insensitive to the screening masses he,m- ^ one strictly assumes r]/s(T = 
T c ) = 0.08, then the ra(T) should be in between and dxj-/dT. On the other hand, 

despite a Ar]/s ~ 0.15 difference in the absolute magnitude, as temperature increases, (i) 
and (ii)’s 77 /s rise at about the same rate, i.e. approximately identical d(r]/s)/dT , and 
both are larger than the one in (iii). Notice that (i) and (ii) have different r^, 7 ie and /im 
but similar /ig — /im near T c . This implies d(r]/s)/dT is sensitive to \ie ~ Hm but is only 
limitedly affected by and the absolute values of /i e,m■ 
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Figure 5. (Color online) (a) The temperature dependence of the dimensionless jet transport 
coefficient q/T 3 for a light quark jet with initial energy E = 10 GeV in the CUJET3.0 framework 
(Red) with scheme: (i) (3.7) (solid), (ii) (3.8) (dashed), and (iii) (3.9) (dotdaslied), compared with 
corresponding q' (Green) as defined in Eq. (3.17). The CUJET2.0 qr/T 3 with {a ma x, Ie, /m) = 
(0.39,1,0) (Blue) and the Af = 4 SYM <isym/T 3 = 71 , 5 ^ V A [145] with ’t Hooft coupling 
A G [ 7 t, 47 t] (Black shaded) are plotted as references, (b) The shear viscosity to entropy density 
ratio 77 /s estimated in the kinetic theory using extrapolation Eq. (3.14) from < 7 ’s in (a). Note 
that there is a clear 77 /s minimum at T ~ 210 MeV in the CUJET3.0 framework regardless of 
the schemes been chosen. The corresponding ( 77 /s)' (determined from q') converges to the pQCD 
weakly-coupled QGP limit at high temperature in (i)(ii) (iii) as expected. The (■ r]/s) m i n in the fast 
liberations always sit above the quantum bound while in the Polyakov liberation it does not. In the 
near T c regime within the fast liberation schemes, the relative magnitude of 77 /s’s does not follow 
the naive inverse of the quark </f’s. This is because the computation of the (ti/s)' receives enhanced 
contributions from softer scales that have stronger electric couplings, and consequently suppressing 
the transverse mean free path. 


3.2.1 Alternative determination of jet transport coefficient 


The q computation above has followed the previous CUJET3.0 prescription [71] as given 
in Eq. (3.10), where the scattering kernel for the jet transport coefficient q is symmetric 
under inter-exchange of E and M in accord with E-M duality considerations. There is 
however a subtle ambiguity: the form of scattering potential in Eq. (3.10) differs from 
the scattering potential in the generalized kernel Eq. (2.15) (as given in the second line of 
Eq. (2.23)) that is used in the actual CUJET3.0 modeling. The l/q|_ factors, present in 
Eq. (2.23) while absent in Eq. (3.10), increase the weight of soft momentum transfers in 
the computation of q. It is important to examine the results for q and rj/s determined from 
the following alternative q' measure of the CUJET3.0 quenching field, and compare them 
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with the computation from Eq. (3.10): 


*< E ’ T, S ° ET <!C,i (qi + /iK 2 W)(qi + /|,f* 2 W)'’ (T) 

x | [C’qqfq + Cqgfg\ ' [«s(ql)] ‘ [/IqI + /1/m^ 2 ( z )] + 

[C 9m (l - f q - f g )] • [1] ■ [/jkqi + /1 /m^ 2 ( z )] } • 


(3.17) 


The “prime” generalizations of quenching parameters of gluon and monopole jets follow 
straightforwardly. By substituting Eq. (3.17) into Eq. (3.14), one can compute the corre¬ 
sponding {p/s)' in the quasi-particle picture according to kinetic theory. Fig. 5(a) shows 
the temperature dependence of both q/T 3 and q' /T 3 for a light quark jet with initial en¬ 
ergy E = 10 GeV in the CUJET3.0 framework with scheme: (i) (3.7), (ii) (3.8), and (iii) 
(3.9). Fig. 5(b) shows the corresponding comparison of 77 /s and {p/s)'. There is a clear 
77 /s minimum at T ~ 210 MeV in the CUJET3.0 framework in both ways of determining 
the quenching parameter. The (77 /s ) 7 nicely converges to the weakly-coupled HTL QGP 
limit at very high temperature T > 500 MeV in all (i) (ii) (iii) schemes, as expected from 
Eq. (2.23) in the \T —> 1 limit. Interestingly, for both estimates of 77 /s, the ( r}/s) m i n in 
the “fast” quark liberation schemes stay above the quantum bound while in the “slow” 
quark liberation scheme it does not. The general relations between [liberation schemes 
+ screening masses] and [(r 7 /s) m , n + d{p/s)/dT} that one could infer from Fig. 4 do not 
alter significantly for the ( 77 /s)' results. Within the fast liberation schemes, in the near T c 
regime, the relative magnitude of 77 /s’s do not follow the naive inverse of the quark qp' s. 
This is understandable since the computation of ( 77 /s )' from Eq. (3.17) puts more weights 
on softer scales that have stronger electric couplings, given ctpctM = 1 at all scales. Thus 
the important EM scattering channel is not affected while the transverse mean free path 
will be suppressed due to larger EE scattering channel cross sections. Consequently one 
gets smaller ( 77 / s)' values as compared with the 77 /s values. 

Overall, the lesson that one learns from Fig. 4 is that rd and he — Hm determines 
(77 / s') min and d{p/s)/dT respectively. Combined with the observations that one draws 
from Fig. 3, we can arrive at the following: 

(1) If data of high pp light hadron (LH) and open heavy flavor (HF)’s Raa and v 2 at 
RHIC and LHC with sufficiently small uncertainties become available, then in the CU- 
JET3.0 framework, after one has constrained the model parameters with LH’s Raa and 
lattice calculations, one can use (1) LH’s V 2 to estimate He{T) — Hm{T ) near T c ; (2) HF’s 
Raa to constrain r^(T); (3) HF’s V 2 to limit the absolute values of he and 7 im as well 
as to further constrain rd(T). Even if the (rd(T), he{T), hm{T)) are not completely fixed 
in CUJET3.0, insights on q(T ) and rj/s(T) can be drawn within this model because for 
jet energy loss in sQGMP, (a) q(T) is regulated by r^(T); (b) d(rj/s)/dT is shaped by 
He(T) - Hm(T) near T c ; (c) {r)/s) min is bounded by r d (T). 

(2) In addition, the investigation of how the near T c physics affects the temperature depen¬ 
dence of the bulk viscosity £/s(T) as well as the role that £/ s(T ) plays on the experimental 
observables at RHIC and LHC is a topic of significant interests, there have been studies 
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from the soft hydro sector [153]. Exploring the importance of (/s(T ) from the hard jet 
quenching sector within CUJET3.0 will be explored elsewhere. 

3.3 Theoretical uncertainties related with monopole density constraints 

Thus far we have concentrated on using the total lattice QCD pressure, p(T), to constrain 
the chromo-magnetic-monopole (cmm) density assuming an ideal gas of chromo-electric- 
charged (cec) and cmm degrees of freedom that leads to 

Pm(T) = i P p{T)/T - p q (T) - p g {T) = p£ S \T) . (3.18) 

Where £ p = 1.012 for a N c = 3, Nf = 2.5 Stefan-Boltzmann gas as in Eq. (3.11). We refer 
to Eq. (3.18) as the Pressure Scheme (PS) to fix the partial pressure of magnetic monopoles 
from the total pressure minus the suppressed semi-QGP densities. 

Such Pressure Scheme may “suffer” from the potential or “bag” B(T ) contribution to 
thermodynamics whereby pressure p = Ts/A — B and energy density e = 3Ts/4 + E>. In this 
case one would then have nonzero “trace anomaly” e—3 p = 4 B and indeed lattice QCD data 
have shown the existence of such a “bag” contribution. In this regard, the entropy density 
s = (e +p)/T, free from any “bag” terms, may serve as a useful “counting” scheme for 
quasiparticle densities. We therefore also introduce an independent Entropy Scheme (ES) 
for determining the total number density via p{T) = £ s s(T)/ 4, where £ s = £ p /4 = 0.253 
for a N c = S,Nf = 2.5 Stefan-Boltzmann gas, and examine the corresponding uncertainty 
in our modeling. 

As can be seen from Fig. 6 , in the Entropy Scheme (ES) the total quasiparticle density 
is higher than in the Pressure Scheme (PS) due to the bag constant: 

p^ s \T) = d s s(T) - Pq {T) - p g (T) = ti s \T) + £ P B(T)/T . (3.19) 

Choosing the ES vs the PS scheme to fix the monopole density would increase the monopole 
density near T c by a factor of 2 and increases the q by approximately 50% near T c . To fit 
the same reference path averaged Raa we would need to adjust the (a c , c m ) in CUJET3.0 
and the (a max , Ie, /m) in CUJET2.0 somewhat respectively. Fig . 8 shows the absolute q 
and the dimensionless q/T 3 in the two schemes. After readjusted (a c , c m ) in CUJET3.0 and 
{otmaxi / e , / m ) in CUJET2.0 to (0.6, 0.33) (as shown in Fig. 7) and (0.35,1, 0) respectively, 
the q near T c in ES is around 50% larger than in the PS. 

The 77 /s in CUJET2.0 and 3.0 computed from inversing the q/T 3 according to Eq. (3.14) 
are plotted in Fig. 9. One sees that the 77 /s in the ES scheme hardly deviates from its value 
in the PR scheme. This is understood since 77 /s is dominated by the “free” quasi-quark 
degrees of freedom. The monopole fluid is almost viscous free in either scheme as one has 
already seen in Fig. 4. 

From this consideration we see that requiring the consistency of hard and soft probes 
can only determine a lower bound on the monopole density near T c . Our default PR scheme 
is above that lower bound. We leave the search for that lower bound to a future study. 
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Figure 6. (Color online) (a) The effective ideal quasiparticle density, p/T 3 = £ p P/T 4 , in the 
Pressure Scheme (PS, Blue) is compared with effective density, p/T 3 = (; p S/4T 3 , in the Entropy 
Scheme (ES, Red) based on fits to lattice data from HotQCD Collaboration [56]. The difference is 
due to an interaction ’bag” pressure —B(T)/T 4 (Green) that encodes the QCD conformal anomaly 
e—3 p ^ 0. (b) The density fraction of the electric (E, red) and magnetic (M, blue) degrees of freedom 
in the \t (solid, Eq. (2T6)) and Xt (dashed, Eq. (3.6)) liberation scheme. The dimensionless E 
and M density p/T 3 in the two schemes are shown in (c )xt an< ^ (d)XT respectively, where both the 
p/T 3 in the PS (solid) and ES (dashed) are plotted. In both liberation schemes, the p m in the ES 
near T c is around twice the p m in the PS. 


4 Conclusions 

We have presented a detailed study of the jet energy loss in semi-Quark-Gluon-Monopole- 
Plasmas (sQGMP), within the new CUJET3.0 framework of jet quenching in bulk con¬ 
strained (VISH2+1D) viscous hydrodynamic backgrounds by extending the perturbative 
QCD based (CUJET2.0) model to include possible non-perturbative chromodynamical fea¬ 
tures of the QCD confinement phase transition near T c ~ 160 MeV. We test the robustness 
and consistency of this new model by comparing predictions of the leading hadron nu¬ 
clear modification factor, Raa(pt > lOGeV/c, y/s), and its azimuthal elliptic asymmetry 
V 2 (pt > lOGeV/c, y/s) with available data from nuclear collisions at RHIC(i/s = 0.2 
ATeV) and LHC(y / s=2.76 ATeV). The sQGMP model depends on two parameters: (1) 
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Figure 7. (Color online) In the CUJET3.0 model with Entropy Scheme (ES) and Xt liberation, the 
a c and c m is adjusted to 0.6 and 0.33 to fit to the reference datum at LHC R’aa(pt = 12.5GeV) « 0.3 
as well as the lattice /xm ([58], c.f. Fig. 2(b)). With this parameter setup, the 7 T°/h ±, s high p p 
Raa and at RHIC Au+Au 200GeV and LHC Pb+Pb 2.76TeV 20-30% collisions are in perfect 
agreements with data [1—9]. The result of prompt D meson, B meson, and heavy flavor decay e— 
results in the ES scheme is plotted in red, green, and orange, respectively. 


the value of the QCD running coupling a c ~ 0.95 — 1.33 at low Q < T c and (2) the ratio 
c m = g^E/R m of nonperturbative electric to magnetic screening scales. We study three 
specific cases, c m = 0,0.3, 0.4. The value of a c is fixed for each case by fitting a single 
reference datum, R c p bPb (pT = 12GeV/c) « 0.3 at LHC at 20-30% centrality. Consistency 
with all available data is then tested comparing predictions of CUJET3.0 on R AA (p p ) and 
v 2 (Pt): f° r h = 7T, D, B at both RHIC and LHC. The emergent chromomagnetic degrees of 
freedom in the sQGMP model near T c are shown to solve efficiently the long standing Raa 
vs V 2 puzzle by leading to a broad maximum of the jet quenching parameter q(E,T) /T 3 
between (1 — 2 )T C . In addition and most remarkably, by extrapolation of the sQGMP q 
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Figure 8. (Color online), (a) The q{T) for a hard quark with initial energy E = 2 GeV 
(dotdashed), 10 GeV (solid), 50 GeV (dashed) computed according to Eq. (3.10) in the Pressure 
scheme, for CUJET3.0 (red), CUJET2.0 (blue) and Af = 4 SYM (black). The dimensionless 
q(T)/T 3 is plotted in (c). (b)(d) The counterpart of (a)(c) in the Entropy scheme. Note that 
( a c ,c m ) in CUJET3.0 and (a max , /e, /m) in CUJET2.0 has been readjusted to (0.6,0.33) and 
(0.35,1,0) respectively fit to the LHC reference datum (cf. Fig. 3). The q in the ES near T c is 
~ 50% larger than in the PS due to the “bag” contribution. 


down to thermal energy E ~ 3T scales one finds the shear viscosity to entropy density 
ratio r)/s ~ T' 3 /q ~ 1/47T to be near the unitary bound, in the critical (1 — 2 )T C transition 
temperature range, which is consistent with viscous hydrodynamic fits to bulk azimuthal 
harmonics v n at px < 2GeV/c. A main goal of the present study is to test the robustness 
of such phenomenological success in CUJET3.0 with respect to a number of theoretical 
uncertainties in the model implementation. 

A key theoretical uncertainty of this model is the rate at which color electric degrees 
of freedom are liberated near the deconfinement transition temperature T c . We have in¬ 
vestigated three very different schemes of color liberation and found that such rate does 
not alter the model’s agreements with high px single light hadrons (LH)’ Raa■ However, 
this rate (r^ = dxx/dT ) considerably influences the temperature dependence of the electric 
screening mass. It is found that the light hadrons’ V 2 (px > 10GeV) is an observable that 
is sensitive to the difference between the electric (E) and magnetic (M) screening mass 
near T c . For open heavy flavors (HF), rd regulates their high px Raa’s hr CU- 
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Figure 9. (Color online), (a) The 77 /s in CUJET2.0 for quasi-quarks (q, red), quasi-gluons 
(g, green), and its total value (black) computed from inversing the q/T 3 according to Eq. (3.14). 
The solid lines correspond to the PS scheme, while the dashed lines correspond to the ES scheme, 
(b) The counterpart of (a) in CUJET3.0. Note that the addition of the monopole (m, blue) like 
quasi-particle degrees of freedom in sQGMP does not alter the overall 77 /s significantly since the 
strong magnetic coupling shrinks the transverse mean free path for monopoles and suppresses the 
shear viscosity contributions from monopoles. Since the sQGMP is dominated by monopole degrees 
of freedom near T c , the total 77 /s in the PS and ES scheme then should naturally converge to the 
same value. 


JET3.0, and their V 2 , s are affected by both the screening masses and the r^. In terms of 
the jet quenching parameter q, it is influenced by rd, but is insensitive to the screening 
masses. On the other hand, for the shear viscosity to entropy density ratio 77 /s, (r]/s) m i n 
and d(r]/s)/dT would vary strongly if rd and [Le — \im changes. 

Therefore, for jet quenching in sQGMP from the CUJET3.0 framework, after model 
parameters are constrained by data of LH’s Raa and relevant lattice calculations, the rate 
of deconfinement rd and the screening masses ^e,m affect jet fragments observables in 
different ways: (1) fiE ~ I^M near T c influences LH’s v% (2) rd influences HF’s Raa', (3) 
rd and absolute values of he,m influence HF’s V 2 - On the other hand, the CUJET3.0’s jet 
transport coefficient q and the shear viscosity 77 /s are affected by and he,m differently: 
(a) q(T) is constrained by r^T); (b) d(r}/s)/dT is constrained by 77 e(T ) — hm(T) near 
T c ', (c) ( 7 ?/s) m i n is constrained by r^T). Given our findings of these dependences, it is 
expected that such model uncertainties can be significantly reduced by experimental input 
such as future high precision data for LH and HF’s Raa as well as vi at both RHIC and 
LHC. For example, the HF’s Raa limits rd(T) therefore constrains q(T) and ( 77 /s) m i n , LH’s 
V 2 limits [1e ~ Hm near T c therefore constrains d(r)/s)/dT, while the two can be combined 
to infer the functional shape of r//s(T); in addition, HF’s V 2 can be utilized to constrain 
the absolute magnitude of fj, e(T ) and 7 : jlm{T). 

Additional discussions on three aspects of the present model have also been included 
in the appendices. The first is about different relativistic corrections from viscous hydro- 
dynamical flows which are found not to affect either Raa or V 2 for light hadrons. This is as 
expected since under eikonal approximation, the number of parton-medium scatterings is 
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fixed in any frame once the initial jet production coordinate and the azimuthal jet propaga¬ 
tion angle are specified. The second aims to identify which one of the various ingredients in 
CUJET3.0 makes the most significant contribution toward obtaining a successful descrip¬ 
tion of the high px azimuthal elliptical harmonics data. We find that the emergent chrorno 
magnetic monopoles play the decisive role. As long as monopoles are present in the near T c 
regime, V 2 is insensitive to the detailed form of the running coupling a s (Q 2 ), provided that 
lattice screening masses can be reasonably described. The last one is about the path length 
(L) dependence of the light/heavy quark energy loss, which is an informative dynamical 
feature of jet quenching models. It is found that both CUJET3.0 and CUJET2.0 converge 
to be around the pQCD/elastic limit at high temperatures T > 400 MeV. However as T 
drops, starting from T ~ 300 MeV, chromo-magnetic monopoles gradually dominate the 
medium. Consequently we find that the energy loss’s dependence on L starts to deviate 
from pQCD/elastic toward the AdS/CFT-like strong coupling behavior. It is noteworthy 
that the “dead cone” suppression is not altered by the nonperturbative sQGMP near T c . 
It is also found that the modification on the L dependence of light and heavy quark energy 
loss caused by the transition from high-T QGP to near-T c sQGMP are quite similar. 

In summary, with our detailed investigations of the jet quenching phenomena in the 
sQGMP model, we conclude that the phenomenological consequences of the sQGMP, for 
both hard and soft probes in heavy ion collisions, stay very robust with respect to certain 
important systematic theoretical uncertainties. The phenomenological success includes 
simultaneous descriptions of all available high px Raa and V 2 data at both RHIC and 
LHC, as well as providing an intuitive dynamical mechanism that may explain how the 
shear viscosity to entropy density ratio could approach the 1/47T unitary bound in the 
vicinity of T c as required to explain the bulk low px < 2 GeV “perfect fluidity”. The 
sQGMP model therefore provides a first model that consistently accounts for both hard 
and soft transport properties of the new state of QCD matter discovered in high energy 
nuclear collisions at RHIC and LHC. 

We end by emphasizing again the fact that remarkably different q(T ) dependence could 
be consistent with the same Raa data, demonstrates clearly the inadequacy of focusing on 
the jet path averaged quantity (q) as the only relevant medium property to characterize 
jet energy loss. Evidently while the {q) captures the important transverse “kick” factor, 
there are other essential factors like the actual chrorno electric and magnetic composition 
of the plasma, the screening masses and the running couplings at multiple scales which 
all strongly influence jet energy loss and imprint their effects beyond just in the (q). It 
is of significant interest and importance to thoroughly exploit these multiple facets of jet 
physics and the opportunities they offer for progressing toward an ultimate understanding 
of the microscopic making of the sQGP. 
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A Relativistic flow corrections to jet energy loss 


Generally speaking, in pQCD based energy loss models, the non-Abelian bremsstrahlung 
amplitudes are derived assuming a static QGP medium [16-20]. With a hydro dynamically 
evolving background, flows move at relativistic velocities in certain cells at certain time; if 
a light-like high-p^ parton happens to pass through a unit cell flowing close to the speed 
of light, since the radiative formulas are applicable only in the co-moving frame, then one 
should boost back to the lab framework for correct predictions of jet quenching observables. 
Therefore, modifications on a pQCD energy loss theory are necessary if it is coupled to a 
hydrodynamically expanding medium. 

In [97, 98], the relativistic flow corrections to q are calculated using weak and strong 
coupling approaches. The authors in those papers derived that in existence of hydrody- 
namical flows, the q should be multiplied by a factor 


r(f) 


Po 


7 f(z) 


1 - 0j(z) • p f (z) 


(A.l) 


Where z and p^ = (po,p) is the coordinate and the four momentum of the high-p-r parton 
in the lab frame, respectively. Note that /3j = p/po- The u j = 7/(1,/?/) is the flow four 
velocity. 

On the other hand, from naive considerations of the jet opacity and relativistic boosts, 
one can get 
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Therefore, to systematically study the relativistic corrections to the energy loss kernel 
hence jet quenching observables, we compare three schemes: (1) T = 7/(1 — fij ■ /3/); (2) 

r = i - /V Pr, (3) r = i. 

Fig. 10(a) (b) shows the results of high pt charged particles’ Raa(pt ) and V 2 {pt) at 
LHC Pb+Pb y/sNN = 2.76 TeV 20-30% centrality collisions using scheme (1)(2)(3) within 
CUJET3.0 and CUJET2.0, compared with available data [5-9]. While all the six combi¬ 
nations generate good agreements with the Raa > only the CUJET3.0 results are consis¬ 
tent with the v r 2 . It is because the nonperturbative sQGMP in the CUJET3.0 framework 
strongly enhances the strength of parton-medium interaction near T c , effectively increases 
late time jet energy loss, and boosts the underestimated V 2 to be in line with data [71]. 
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Figure 10. (Color online) Charged particles’ (a) Raa and (b) i >2 in LHC Pb+Pb yijv\iv=2.76TeV 
semi-peripheral collisions, computed from CUJET2.0 HTL fs = 1, /m = 0 (purple) and CUJET3.0 
c m = 0.3 (red) with relativistic flow corrections (1) T = 7/(1 — /3 ] ■ /3/) (solid) [97, 98]; (2) P = 
1 — 0j ■ Bf (dashed); (3) P = 1 (dotdashed) to the energy loss kernel, compared with relevant 
data [5-9] (blue). The parameters a max (v2.0) and a c (v3.0) are adjusted to fit to the R aa (pt = 
12.5GeV) « 0.3 reference datum. Note that both Raa and Vi do not distinguish the different 
flow corrections at a measurable level. This can be partially understood as the number of parton- 
medium scatterings is fixed for a given jet path in any frame once the initial production coordinate 
and azimuthal propagation angle have been specified. 


Surprisingly, as shown in Fig. 10(b), the relativistic flow correction scheme (1)(2)(3) 
do not alter the V 2 significantly in both CUJET3.0 and CUJET2.0, despite very different 
values of a c and a max (note that a c ’s and a mai ’s are adjusted to fit reasonably the reference 
datum R^ a {pt = 12.5GeV) ~ 0.3). Under eikonal approximation, a jet path is fixed once 
the initial jet production coordinate and the azimuthal jet propagation angle are specified. 
Therefore, the robustness of the Raa and V 2 in response to relativistic flow corrections 
can be understood as along the fixed jet path, the number of parton-medium scatterings 
is invariant in any frame. To be more careful, one would argue that the V 2 shifts by ~ 0.01 
from scheme (3) T = 1 to (1) T = 7/(1 — (3/ • /3/), which is not “unchanged”. In fact, this 
magnitude of variation is consistent with the conclusion of [98], where the authors argue 
that the corrections on ordinary jet quenching observables because of the hydro flow are 
too small to be measured hence can be neglected at current stage. 

B What contributes most significantly to the strongly enhanced jet opac¬ 
ity in sQGMP? 

In CUJET3.0, the dimensionless jet quenching parameter g/T 3 is strongly enhanced in 
the near T c regime, and several factors may contribute to this enhancement: the enlarge¬ 
ment of a s —> a c as T —> T+ in Eq. (2.24), the separation of the electric and magnetic 
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Figure 11. (Color online) (a) The running strong coupling a s (Q 2 ) in four different models, 
[a c /a maa ,]-t-[QGMP/HTL]. Notice that a max (dashed) and a c (solid) is parametrized as in Eq. (2.3) 
and Eq. (2.24) respectively. Note in CUJET3.0 energy loss kernel, i.e. Eq. (2.23) [QGMP], both 
chromo-electric (red) and chromo-magnetic (blue) coupling exists, while in CUJET2.0 energy loss 
kernel, i.e. Eq. (2.1) [HTL], only the electric coupling (green) is present, (b) The electric screening 
mass (he) and magnetic screening mass (pm) as temperature varies in the four models, compared 
with lattice data (E, red; M, blue; [58]). Note in plotting the [u c /a max \ + [QGMP/HTL] results, 
the curve styles are the same as in panel (a). In [a c ], T c = 160 MeV; In [a mox ], A qcd = 200 MeV. 
Note that in [HTL] the magnetic screening mass is zero because /m = 0 [95]. 


quasi-particle density fraction according to Eq. (2.16), and the magnetic screening mass 
regulator c m in Eq. (2.21). A critical question to ask is: which factor contributes most 
significantly to the enhancement of the jet opacity in sQGMP? In order to answer this, 
a practical solution is to hybridize the CUJET2.0 running coupling scheme (denote it as 
Oimax-, Eq- (2.3)) with the CUJET3.0 energy loss kernel (denoted it as QGMP, Eq. (2.23)); 
and hybridize the CUJET3.0 running coupling scheme (denote it as a c , Eq. (2.24)) with 
the CUJET2.0 energy loss kernel (denot it as HTL, Eq. (2.1)); then compare the predic¬ 
tions of jet quenching observables, in particular, high px Raa and v 2 from the four models: 
[etc/Umax] + [QGMP/HTL], 

Fig. 11(a) shows the running coupling used in the four models, note that the parameters 
are fixed by fitting the LHC R!^ a {pt = 12.5GeV) « 0.3 reference datum, c.f. Fig. 12(a). 
Note that the [a maa; ]-|-[QGMP] model has a weaker chromo-electric coupling in the near 
T c regime but a stronger one at large Q > 2 GeV than the [a c ]+[QGMP] model. This 
model also has a smaller otf nax = 0.24 than the a ma x = 0.39 in the [a: c ] + [HTL] model. It is 
as expected since [a max \ +[QGMP] has an extra monopole fraction with extremely strong 
couplings. 

Fig. 11(b) shows the electric and magnetic screening masses in the four models com¬ 
pared with lattice data [58]. Both [a max ]+[QGMP] and [ct c ] + [QGMP] can describe both 
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Figure 12. (Color online) Charged particles’ (a) Raa and (b) i >2 in LHC Pb+Pb y / sjviv=2.76TeV 
semi-peripheral collisions, computed from CUJET2.0 [HTL] /e = 1, /m = 0 (green, Eq. (2.1)) and 
CUJET3.0 [QGMP] c m = 0.3 (red, Eq. (2.23)) with the [ a max \ (dashed, Eq. (2.3)) and [ay] (solid, 
Eq. (2.24)) running coupling scheme, compared with available data [5-9] (blue). The parameter 
oimax and a c are adjusted to fit to the R\ a {Pt = 12.5 GeV/c) ~ 0.3 reference point. Note that 
while all four models can explain the high pt Raa , only [a c ] + [QGMP] and [a m0 i] + [QGMP] can 
explain the high px V 2 - This suggests the emergence of chromo-magnetic monopoles as T —> T c + 
contributes most significantly to the strongly enhanced q/T 3 near T c and generates the simultaneous 
description of high px light hadrons’ Raa and V 2 ■ 


Pe and pm reasonably well. Both [a maa; ]+[HTL] and [a c ]+[HTL] are in agreements with 
Pe, but they have pm = 0 because /m = 0 in [HTL] [95]. 

Fig. 12 show the results of high px light hadrons’ (a) Raa and (b) V 2 in the four models, 
i.e. [ay or a max \ + [QGMP or HTL], at LHC Pb+Pb ^snn = 2.76 TeV 20-30% collisions, 
compared with available data [5-9]. Note the model parameters in the four models are 
fixed as: 


[ay] + [QGMP] : a c = 0.95, c m = 0.3 ; (B.l) 

[a max ] + [QGMP] : a max = 0.24, c m = 0.3 ; (B.2) 

[a c ] + [HTL] : a c = 2.38, f E = TO, f M = 0.0 ; (B.3) 

[a max ] + [HTL] : a max = 0.39, f E = TO, f M = 0.0 . (B.4) 


While all four models are compatible with the Raa data, only [ay or a max ] + [QGMP] can 
fit to the high px charged particles’ V 2 - Noted that the boost in azimuthal elliptical har¬ 
monics is contributed mainly by the enhancement of the jet opacity near T c ; [a c ]+[QGMP] 
and [a maa; ] + [QGMP] have different running coupling schemes as well as different pm 's in 
the near T c regime but share the same CUJET3.0 energy loss kernel; we can therefore con¬ 
clude that the dividing of electric (E) and magnetic (M) fractions according to Eq. (2.16) 
results in the transition of the nonperturbative medium from E dominate to M dominate as 
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temperature goes towards Tj 1 ", and contributes most significantly to the strongly enhanced 


parton-nredium interaction near T c for the jet energy loss in semi-Quark-Gluon-Monopole- 
Plasmas. 

C Path length dependence of parton energy loss in sQGMP 

CUJET3.0 is a jet quenching framework that hybridizes the perturbative dynamical DGLV 
opacity expansion theory, the TG elastic energy loss, and the nonperturbative sQGMP near 
T c . It is consistent with high px Raa and V 2 data at RHC and LHC simultaneously, and 
intrinsically builds a connection between high energy jet quenching and bulk perfect fluidity 
[71]. Beyond this phenomenological success, a crucial question to ask is how does the path 
length dependence of jet energy loss look like in this hybrid model. From addressing this, 
one can gain insights into e.g. at what temperature does the nonperturbative physics enter 
jet quenching, why different q(T )’s lead to the same suppression factor in CUJET3.0 and 
CUJET2.0, whether or not the beyond leading order effects change light and heavy quark 
energy loss identically, etc. 

To obtain useful insights about the path length dependence of jet energy loss in general 
scenarios, one can take a step back to look at the so-called “abc” model [27, 28] where the 
parton energy loss is simplified as a power law of the energy E, the path length L, and the 
local temperature T: 



i a t brp2—a-\-b 


(C.l) 


Depending on underlying energy loss mechanism, the b may take quite different values 
(see discussions in e.g. [27, 28]: for collisional processes dominated energy loss, b m 0; at 
leading order (LO) in pQCD, for non-Abelian bremsstrahlung processes dominated energy 
loss, b « 1; as the coupling strength a s becomes extremely strong (as well as N c -A 
oo) and AdS/CFT correspondence is applicable, holographic jet quenching generally has 
b ~ 2. In the following subsections, after we have interpolated (A E/E)(L) at various 
temperatures, we will further extract out a “b” factor from b = dlog(AE / E) / dlog(L / Lq) — 
1. Different from the abc model with “global” power law dependence for the jet energy 
loss, our extracted “b” factor would be a sort of “local index” but it nevertheless is an 
informative indicator that can help achieve deeper understandings about the energy loss 
dynamics encoded in the computed (A E/E){L) from the CUJET3.0. 

C.l Light Quark 

In the radiative sector, there have been next-to-leading order (NLO) pQCD calculations for 
energy loss assuming massless projectile partons [154-157], and they all suggest a double 
logarithmic path length dependence of the jet quenching parameter, i.e. 


q oc log L( 1 + log 2 L) . 


(C.2) 


Where L = L/Lq and Lq is a proper ultraviolet cutoff. Generally speaking, the differentiate 
jet energy loss has dE/dL oc qL. Let us assume logL varies much slower than L, after 
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straightforward integrations and simplifications, one arrives at 

log(A E/E) oc 2 logZ + log log L + log(l + log 2 L) + const . (C.3) 

This form is of cause not general enough because of the logL approximation and it is 
derived at NLO in pQCD. Nevertheless, motivated by Eq. (C.3), we will use 

log —^ = A + B log + log(l + C log ^-) + log(l + D log 2 ) , (C.4) 

to interpolate the path length dependence of radiative jet energy loss in sQGMP within 
the CUJET3.0 framework. On the other hand, combining the above with Eq. (C.l), one 
can easily get 

= dlog(AEy E) 
d\og{L/Lo) 

ft _ 1 I __I_^_ 

l + Clog(L/L 0 ) l + Dlog 2 (L/L 0 ) 

—> B — 1 (L —> oo) . 

The upper panels of Fig. 13 show the path length dependence of the radiative energy 
loss A E/E of a high-p^ light quark (mass M = 200 MeV) with initial energy E = 20 
GeV traversing a brick plasma (fixed density) at varied temperatures in CUJET3.0 and 
in CUJET2.0. The (a)(b)(c)(d) corresponds to jet quenching in the brick medium with 
temperature T = 450, 300, 200, 160 MeV respectively. In the computation, Eq. (2.23) 
and Eq. (2.1) is used for CUJET3.0 and CUJET2.0 respectively, and the Monte-Carlo 
integration is iterated 1,000,000 times to enforce convergence. The brick size L is increased 
from 0.45 fm to 3.05 fm with 0.5 fm intervals. We fit the A E/E vs L in both CUJET3.0 
and CUJET2.0 with Eq. (C.4). Note the phenomenon that A E/E exceeds 1 for some large 
size bricks is purely technical, since in the computation E is not dynamically updated in 
Eq. (2.23) and (2.1) as hard partons propagate through and lose energies. Nevertheless, 
the A E/E(L) at small L’s, the slope of AE/E(L), the thickness where A E/E = 1, and 
the relative information about parton energy loss at the four different temperatures are all 
meaningful. 

For the CUJET2.0 model with pQCD+HTL, as temperature decreases, d(AE/E)/dL 
is altered significantly. Define the stopping distance L\ as 

Li : AE/E\ L=Ll = 1 . (C.8) 

It monotonically increases in this picture. This is as expected, since in the CUJET2.0, 
bricks with lower temperatures have less opacities, and it takes a longer path for a jet to 
lose all the energy. However, for the L\ in CUJET3.0, though it monotonically increases, 
its increasing rate is much less than CUJET2.0’s. More importantly, the CUJET3.0’s 
d(AE/E)/dL significantly grows as the temperature decreases below T = 300 MeV. 
All these pointed to the fact that, as the temperature approaches T c , chrorno-magnetic 
monopoles gradually dominate the medium, since olm oce, the jet-medium interaction 


(C.5) 

(C.6) 

(C.7) 
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Figure 13. (Color online) (Upper panels) The radiative energy loss ratio A E/E of a light 
quark jet (M = 200 MeV) with initial energy E = 20 GeV propagating through a brick plasma 
with various thicknesses L at temperature T = (a) 450, (b) 300, (c) 200, (d) 160 MeV, in the 
CUJET3.0 a c = 0.95, c m = 0.3, Xt model (red) and in the CUJET2.0 a max = 0.39, /# = 1, 
/e = 0 model (blue). When interpolating the AE/E(L)’s , the NLO pQCD motivated Eq. (C.4) 
is used. Notice that A E/E exceeds 1 at some large L for a certain T, this is due to the fact that 
in the computation E is not dynamically updated in each unit cell. Since with a smooth viscous 
hydro evolution background, the hard parton will stay at a certain temperature for limited time, 
the relevant energy loss informations are mostly at short L. As temperature decreases, CUJET3.0’s 
A E/E(L) and the stopping distance L\ (defined in Eq. (C.8)) respectively gets steeper and larger 
compared with CUJET2.0’s. At low and intermediate T, (A£/£)„ 3 .o < (A E/E) V 2 .o in the small 
L regime. (Lower panels) The path length L dependence of the power b in Eq. (C.5) at different 
temperatures. Note that b = 0, 1, 2 is approximately the elastic, pQCD and AdS limit respectively. 
At high temperature T ~ 400 MeV, the b(L) of CUJET3.0 and CUJET2.0 converge to around the 
pQCD limit. As temperature cools down, when T ~ 300 MeV, CUJET3.0’s b(L) start becoming 
larger than CUJET2.0’s. This signals the transition from E to M dominant as well as from weak to 
strong coupling for the bulk. In the near T c regime, the b(L) in the CUJET3.0 framework is higher 
than LO pQCD, and is close to the AdS limit. This implies the model ingredients in CUJET3.0 do 
effectively bring in nonperturbative dynamics into the original pQCD/DGLV energy loss kernel. 


may weaken at a reduced rate or it may be enhanced (telling from Fig. 4(a), this is the 
correct picture) despite the decrease of the quasi-particle number density. 

To this point, there is a critical question that one should address regarding the jet en¬ 
ergy loss in sQGMP: if one compares the CUJET3.0 and CUJET2.0’s q/T 3 in Fig. 4(a)(b), 
the former is always above the latter in the temperature range of T < 450 MeV, then 
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why both of them can reasonably describe the high pt light hadron’s Raa at RHIC and 
LHC? The upper panels of Fig. 13 give one the answer: at T = 400 MeV, CUJET3.0 and 
CUJET2.0’s A E/E(L) almost overlap; as the temperature cools down, because of the tran¬ 
sition of the nonperturbative medium from EQPs dominate to MQPs dominate, the CU- 
JET3.0’s A E/E(L) becomes steeper and steeper than the CUJET2.0’s, while the former’s 
L\ becomes less and less than the latter’s; these effects lead to (A E/E) v3 _q < (AE / E) V 2 .q 
at small L < L eq , where L eq is defined as 

L eq : (AE/E) v3 , 0 \ L=Leq = ( AE/E) v2 .o\L=L eq , (C.9) 

and this L eq keeps enlarging as T decreases. For hadron suppressions in a hydro dynamically 
evolving smooth medium in A+A collisions, along a given jet path (let the initial production 
point be at the origin), if one sequentially divides it into sections with average temperature 
of 450, 300, 200, 160 MeV and mark the traveling time in each section as Z 450 , Z 300 , ^ 200 > 
Zi60) in a general estimation, the ordering is Z 450 > ^300 > ^200 > Zi60) and Zi60 “C 1 fm while 
Z 200 0.5 fm. Define 

A e(T = -AZ_!±!) = - / dL [{AE) v3 . 0 - (A E) v2 . 0 \ . (C.10) 

Z — L k J r p(L < )=Ti 

Although Ae(T = 450, 300 MeV) are most likely positive, A e(T = 200,160 MeV) can 
be nontrivially negative and thus compensate the formers. This cancellation results in 
the overall A^ « 0, meaning similar averaged Raa predictions for CUJET3.0 and 
CUJET2.0. 

The lower panels of Fig. 13 show the extracted power b (b = d\og{AE / E) / d\og{L / Lq)— 
1, c.f. Eq. (C.5)) versus the brick thickness L. A first observation is that, at high tem¬ 
perature, the CUJET3.0’s b(L) is almost identical to CUJET2.0’s; but as T goes down, 
starting from T ~ 300 MeV, the former gets larger than the latter. This is a clear sig¬ 
nal showing that the chromo-magnetic monopoles begin emerging and bringing up strong 
coupling effects from T ~ 300 MeV ~ 1.8 T c . Let L 2 be the path length that satisfies 

b{L = L 2 ) = 1 . (C.ll) 

In both CUJET2.0 and CUJET3.0, L 2 increases as the temperature decreases, this is 
understood as the opacity L/X = paL ~ T ? 'L , if T drops, a larger L is required to keep the 
same opacity level and hence the similar antenna structure. Meanwhile, the CUJET2.0’s 
db/dL does not undergo significant shifts as temperature varies, but this is not the case for 
CUJET3.0. Define 


61.5 = b(L = 1.5 fm) . (C.12) 

In CUJET2.0, 61.5 rises from ~ 0.6 at 450 MeV to ~ 1.2 at 160 MeV, which is roughly 
consistent with the LO pQCD expectation of b = 1. In CUJET3.0, the 61.5 rises from ~ 0.6 
at 450 MeV to ~ 2.5 at 160 MeV, suggests that the sQGMP introduces some nonperturba¬ 
tive effects into the DGLV opacity series and effectively causes the resummation of higher 
orders in the full QCD amplitude. 
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Figure 14. (Color online) The radiative energy loss ratio A E/E and the power b (c.f. Eq. (C.5)) 
of a high-pT heavy quark (M = 4.75 GeV) with initial energy E = 20 GeV traversing a brick 
plasma with thickness L at various temperatures in CUJET3.0 and in CUJET2.0. All marks and 
computational details are the same as in Fig. 13. Note that the heavy quark’s d(AE / E) / dL and 
b(L) are smaller than the light quark’s as expected from the dead cone suppression. At high T ~ 
450 MeV, both CUJET3.0 and CUJET2.0 converge at around the linear elastic energy loss limit. 
As T drops towards T c , beginning from T ~ 300 MeV, the CUJET3.0’s b(L) starts deviating from 
CUJET2.0’s because of the emergence of chromo-magnetic monopoles; the former’s A E/E(V) gets 
steeper than the latter’s, while L i (c.f. Eq. (C. 8 )) gets shorter. All these alternations for the heavy 
quark suppression are similar to those for the light quark, and the magnitude of the b(L) deviation 
for the two different flavors are almost identical. This suggests that the nonperturbative effects in 
the near-T c sQGMP modify the energy loss kernel of light and heavy quarks in a very similar way. 


C.2 Heavy Flavor 

Let us now turn to the path length dependence of the heavy quark energy loss in the 
sQGMP. The upper panels of Fig. 14 show the A E/E(L) of a heavy quark jet (mass 
M = 4.75 GeV) with initial energy E = 20 GeV transversing a brick plasma at T = 
450, 300, 200, 160 MeV, in CUJET3.0 and in CUJET2.0. Except for the jet mass M, all 
technical details in this computation are the same as the in the one for the light quark 
energy loss. Compared with Fig. 13, one notice that slope of A E/E(L) in both CUJET3.0 
and CUJET2.0 are more gentle for the heavy quark than for the light quark, and the L\ 
(c.f. Eq. (C.8)) grows faster when cooling down. This clearly indicates the dead cone effects 
suppress the induced radiation regardless of whether or not the sQGMP is present. As the 
temperature gets lower, for heavy quarks, CUJET3.0’s A E/E(L) and L\ also becomes 
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steeper and smaller than CUJET2.0’s. This phenomenon has the same physical origin as 
for light quarks discussed in section C.l, i.e. a transition from EQPs to MQPs dominate. 
Interestingly, in the near T c regime, the L eq (c.f. Eq. (C.9)) for heavy quarks is smaller 
than for light quarks. Based on the discussions in section C.l, this will lead to the high pp 
leading B meson Raa predictions from CUJET3.0 being slightly lower than CUJET2.0. In 
fact, this is case from the comparison of the R in [71] and in [38]. 

The lower panels of Fig. 14 show the extracted power 6 (c.f. Eq. (C.5)) versus the 
medium thickness L. Notice that for some temperatures at large L, the b becomes less 
than 1, nevertheless this can be neglected since in these regimes the A E/E has became 
larger than 1, which is unphysical. In CUJET2.0, the 61.5 (c.f. section C.l) rises from 
~ 0.3 at 450 MeV to ~ 0.6 at 160 MeV, this is weaker than the LO pQCD radiative energy 
loss 6 = 1, and approaches the elastic limit 6 = 0. Compared with the energy loss for light 
quarks, the dead cone suppression is significant for heavy quarks. In CUJET3.0, the 61.5 
rises from ~ 0.4 at 450 MeV to ~ 1.6 at 160 MeV, this suggests the strong coupling effects 
hence high order resummations also enter the energy loss kernel for the heavy quark. 

At high T ~ 450 MeV, the CUJET3.0’s and CUJET2.0’s b(L) overlap, as T drops, 
beginning from T ~ 300 MeV, the former starts to deviate from the latter, suggesting 
the commencement of monopoles taking control of the medium transport properties. This 
initiating temperature TimtiaZ ~ 300 MeV for the heavy quark coincides exactly with the 
Tinitial for the light quark; meanwhile, the magnitude of the deviation in b(L) between CU- 
JET3.0 and CUJET2.0 for the two different flavors resemble each other; these observations 
imply that the nonperturbative effects brought up by the sQGMP near T c influences the 
light quark jet quenching and the heavy quark jet quenching in approximately the same 
manner within CUJET3.0. 
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